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Abstract. Quantum mechanical systems exhibit an inherently probabilistic 
nature upon measurement which excludes in principle the singular direct ob- 
servability continual case. Quantum theory of time continuous measurements 
and quantum prediction theory, developed by the author on the basis of an 
independent-increment model for quantum noise and nondemolition causal- 
ity principle in the 80's, solves this problem allowing continual quantum pre- 
dictions and reducing many quantum information problems like problems of 
quantum feedback control to the classical stochastic ones. Using explicit in- 
direct observation models for diffusive and counting measurements we derive 
quantum filtering (prediction) equations to describe the stochastic evolution of 
the open quantum system under the continuous partial observation. Working 
in parallel with classical indeterministic control theory, we show the Markov 
Bellman equations for optimal feedback control of the a posteriori stochas- 
tic quantum states conditioned upon these two kinds of measurements. The 
resulting filtering and Bellman equation for the diffusive observation is then 
applied to the explicitly solvable quantum linear-quadratic-Gaussian (LQG) 
problem which emphasizes many similarities and differences with the corre- 
sponding classical nonlinear filtering and control problems and demonstrates 
microduality between quantum filtering and classical control. 
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1. Introduction 

With technological advances now allowing the possibility of continuous moni- 
toring and rapid manipulations of systems at the quantum level [U [2] , there is an 
increasing awareness of the importance of quantum feedback control in applica- 
tions of quantum information such as the dynamical problems of quantum error 
corrections and quantum computations. The theory of quantum feedback control 
based upon a classical stochastic process formally developed by the author in the 
80's [3j IH [5] has been recently applied in many contexts including state prepara- 
tion [3 [8] , purification [9l [10| , risk-sensitive control [HJ [12] and quantum error 
correction [13l [14] . It has also been studied from the practical point of view of 
stability theory [15] which contains a useful introduction to quantum probability 
and along with ,16] gives a comprehensive discussion on the comparison of classical 
and quantum control techniques. 

The main ingredients of quantum control are essentially the same as in the 
classical case. One controls the system by coupling to an external control field 
which modifies the system in a desirable manner. The desired objectives of the 
control can be encoded into a cost function along with any other stipulations or 
restrictions on the controls such that the minimization of this cost indicates opti- 
mality of the control process. We restrict ourselves to the more interesting case of 
quantum closed-loop or feedback control based on the indirect continuous in time 
observations, the quantum stochastic theory of which was initiated by Belavkin 
in a series of increasing generality papers [TTj, HH [TH [201 as a quantum anal- 
ogy to the classical stochastic prediction theory which is based upon the nonlinear 
(Stratonovich) filtering equations. This work was developed then in the beginning 
of 90 's [2H [221 [231 EH, [25J ESJ [27] , and this development laid down the foundations 
for Eventum Mechanics, a new quantum stochastic mechanics of continual observa- 
tions giving a microscopic theory of continuous reductions [551 HH] and spontaneous 
localizations [30l[3T] for quantum diffusions [32J [33] , quantum jumps [3H[35], and 
other mixed stochastic quantum trajectories [36l [37l [38] . 

In order to demonstrate the power of this new event enhanced quantum me- 
chanics, it was also applied right from the beginning [3l |39l [5] to solve the typical 
problems of quantum feedback control in parallel to the work on classical stochas- 
tic control with partial observations first introduced by Stratonovich 001 01] and 
Mortensen 142!. Thus, the problem of optimal quantum feedback control was sepa- 
rated into quantum filtering which provides optimal estimates of the quantum state 
variables (density operators) and then a classical optimal control problem based on 
the output of the quantum filter. The classical noise which is filtered out by passing 
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from the prior to the posterior quantum states comes from the irreducible distur- 
bance to the quantum system during the observation (due to the interaction with 
measurement apparatus). Unlike in classical case, this is an unavoidable feature 
of quantum measurement since the state of an individual quantum system is not 
directly observable. However, the lack of urgency for such a theory and the novelty 
of the mathematical language at the time left this work relatively undiscovered only 
to be rediscovered recently in the physics and engineering community. 

The purpose of this paper is to build on the original work of the author and 
present an accessible account of the theory of quantum continual measurements, 
quantum causality and predictions and optimal quantum feedback control. Firstly 
we introduce the necessary concepts and mathematical tools from modern quantum 
theory including quantum probability, continuous causal (non-demolition) mea- 
surements, quantum stochastic calculus and quantum filtering. Next the quantum 
Bellman equations for optimal feedback control with diffusive and counting mea- 
surement schemes are informally derived. The latter results were first stated in 
[5] without derivation and a consideration for the diffusive case was recently given 
in [43]. We conclude with an application of these results to the multi-dimensional 
quantum Linear-Quadratic-Gaussian (LQG) problem and a discussion on the com- 
parison with the corresponding classical results. However we first start from a 
model example of LQG quantum filtering and feedback control problem which is 
important since it is one of the few exactly solvable control problems which em- 
phasizes the similarities between the corresponding classical and quantum filtering 
and control theories. It allows us to set up notations and clearly demonstrates not 
only the similarity but also the difference of classical and quantum feedback control 
theories which can be observed in microduality principle, a more elaborated duality 
between quantum linear Gaussian filtering and classical linear optimal control. 

1.1. Model example: quantum free particle. The quantum linear filtering and 
optimal quadratic control problem with quantum Gaussian noise was first studied 
and resolved by the author in a series of quantum measurement and filtering papers 
[17] , [44] , [28] and based on this quantum feedback control papers [H] , [20] , [5] . The 
simplest example of a single quantum Gaussian oscillator matched with a transmis- 
sion line |44j as complex one-dimensional channel was taken as a quantum feedback 
model in the starting preprint [18] . eventually published in [6]. However a more 
similar to the classical case quantum linear models require at least two real dimen- 
sions instead of single complex, and we may now use the multidimensional quantum 
LQG control solutions derived in the last Section of this paper for application on 
higher dimensional systems which do not have such complex representation. The 
optimal control of a continuously observed quantum free particle with quadratic 
cost is the simplest such example. 

Let x, — {xi,±2) be a pair of phase space operators x\ — q, ±2 — p for a 
quantum particle in one dimension, given by selfadjoint operators of position q and 
momentum p satisfying the canonical commutation relation (CCR) 

(1.1) [q,p] := qp-pq = ifrl. 

Here I is the identity operator in a Hilbert space f) of the CCR representation 
([Lip and h is called Planck constant, which could be for our purpose any positive 
constant h > 0. Let us denote the row of initial expectations (xj) of Xj in a quantum 
Gaussian state by x % = (q, p) , and also denote the initial dispersions of q and p by a q 
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and cr p respectively and the initial symmetric covariance Re {qp) ~ qp by 

® QP — ^ PQ' 

The Hamiltonian p 2 /2/i of free particle is perturbed by a controlling force using 
the linear potential <f>(t, q) = f3u (t) q with u (t) 6l as H (u) = p 2 /2fi + (3uq where 
fi > is the mass of the particle. The particle is assumed to be coupled not only to 
control which can be realized by a quantum coherent (forward) channel, but also to 
a coherent observation (estimation) quantum channel such that its open Heisenberg 
dynamics is described by quantum Langevin equations as a case of f|5 . 3|) : 

(1.2) dQ (t) + XQ (t) dt = -P (t) dt + dW* Q(0) = q 

/i 

(1.3) dP (t) + XP (t) dt = dVp* - (3u (t) dt, P(0)=p 

Here A = \ (ae + pj) and V* = aV* + fiVj, W* = — eW* — "fWj are given by two 
independent pairs (V ,W ) (where o = e, / stands for error and force) of Wiener 
noises V a — Ks(A+), W Q — 2K(A+) due to the interaction with the coupled 
estimation and feedback channels. Note that these noises do not commute, 

(1.4) w;v; - v;w; = (r a s) mi, 

if A 7^ 0, which is necessary and sufficient condition for preservation of the CCR 
(jl.ip by the system (II. 2|) , (|1.3|) . It can be easily found by substituting the solution 

P (t) = e~ xt p + f e (s "* )A (dV p s - (3u (s) ds) 
Jo 

of the second equation (|1.3p into the first (|1.2p that 

[Q(r),Q(s)} = -\r- S \e- x ^- s ^0. 
M 

Therefore the family {Q (t)} is incompatible, cannot be represented as a classical 
stochastic process and directly observed. However it can be indirectly observed 
by continuous measuring of the coupling operator aq with an error white noise 
in the estimation channel as it was suggested in [5].|25|. To this end we measure 
Wl = 25R(A+) as an input process evolved after an interaction with the particle 
into an output classical process given by a commutative family [Y* : t > 0] in the 
linear estimation channel 

(1.5) dy e * = aQ (t) dt + dW*, 

in which the input process appears as measurement error noise with commutative 
independent increments dW* representing the standard Wiener process such that 
(dHy = dt, but noncommuting with the perturbative force V* since, as it is 
explained in the last Section, 



(1.6) dW*dV* = ^-dt, dW*dW* = -edt. 

Thus the measurement error noise W* satisfies the error-perturbation CCR 

(1.7) [Vr,W°]={rAs)ihaI, 

which is necessary and sufficient condition of quantum causality (or quantum non- 
demolition condition) in the form 

(1.8) [Y:,Q(s)}=0=[Y e r ,P(s)} Vr<s 

requiring the statistical predictability of quantum hidden in the future trajecto- 
ries {X, (s) : s > t} with respect to the classical observed in the past trajectories 
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{YJ : r < t} for each t. From this we derive the Heiscnberg error-perturbation un- 
certainty principle in the precise Belavkin inequality form AA\ , 28J 



2 ( cth x 2 



(i.9) wr > (^ T j d*, (d^) z =at 

in terms of the perturbation V p in (|1.3[) and standard error W* in (|1.5[) . Thus we 
have the case 



(1-10) J=l _° x J J, A= I | 



h 

n. 

n 



of the general quantum linear open system considered in the last Section, where 
A is the direct sum A e © Aj of two rows A e , Af corresponding to B e = (a,0), 
E = (0, e), B = (/3,0), Ej = (0,7). From this we compute the matrices (|5.10l) and 
(|5.24p satisfying the microduality principle, which are turned to be diagonal, 



C P J ' v % 

with eigenvalues ( q = 7 2 , ( p = (^/2) 2 (a 2 + fl 2 ) = r) q , i] p = e 2 . 

1.2. Quantum feedback control example. We can now apply the results ob- 
tained in the last Section to demonstrate optimal quantum filtering and optimal 
feedback control and their microduality on this model example. The optimal es- 
timates of the position and momentum based on a nondemolition observation of 
free quantum particle via the continuous measurement of Y t , originally derived in 
5j, 25j in absence of control channel, arc then given by the Belavkin Kalman filter 
(|5.12p in the form of linear stochastic equation 

(1.11) d^ + ^Adi = -pldt + (aa q (t) - e) dW* 

A* 

(1.12) dpl+plXdt = (3u(t)dt + aa qp {t)dWl, 

where the estimation innovation process W\ describes the gain of information due 
to measurement of Y* given by 

(1.13) dWt = dF e * - aq^dt, 
and the error covariances satisfy the Riccati equations 



dt CT <7 



(l-l 4 ) ii^qp = \°p - (A - S) a qp - a 2 a q a qp 

Tt a P = C P - 2Act p - (ao p f, 

where we denote 6 = i (ae — j/3), with initial conditions 

a q (0) = a q , a qp (0) = a qp , a p (0) = a p . 

The Riccati equations for the error covariance in the filtered free particle dy- 
namics have an exact solution [25 with profound implications for the ultimate 
quantum limit satisfying the Heisenberg uncertainty relations for the accuracy of 
optimal quantum state estimation via the continuous indirect quantum particle 
coordinate measurement. 
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The dual optimal control problem can be found by identifying the corresponding 
dual matrices from the table (|5.23p which give the quadratic control parameters 



(1.15) 



c(u) 



s = uj q q 2 + LUgpipq + qp) + u p p 2 

corresponding to the dual output process given by z — "fp. For the linear Gaussian 
system (|5 . 19[) gives the optimal control strategy 

(1.16) u(t)=[3(cj pq (t)p t +u p (t)q t ) 

where the coefficients are the solutions to the Riccati equations 



(1.17) 



(*) 



% - 2Xa q - (f3a q ) 2 

-jW, - (A + S) LO q p - [PtjJpLO 



^1 

r/„ + 2 



+ 2 ^LOqp ~ UpS^J - {flCiJp) 



p^ qp 



with terminal conditions 
uj p (T) 



0J o . 



'pi Uqp {T) — Uqp, (jJq (T) — uj q . 

Note that in this example, as well as identifying the dual matrices by transposition 
and time reversal according to the duality table (|5.23p , one must also symplecticly 
interchange the phase coordinates (q, p) <-> (p, q) . This is because the matrix of 
coefficients A is non-symmetric and nilpotent, so it is dual to its transpose only 
when we interchange the coordinates in the dual picture. Thus the optimal coeffi- 
cients {ujp, LOqp, tOq} (t) in the quadratic cost-to-go correspond to the minimal error 
covariances {<J q , u qpi cr p } (T — t) in the dual picture. 

The minimal total cost for the experiment can be obtained from (|5.22p by sub- 
stitution of these solutions 



S = aj q (q 2 + <j q ) + 2ujqp(qp + a qp ) 

+u p _(Q) ( P 2 + op) + ^{h 2 uj p (t) + to 2 pq (t) o q (t))dt 
2u n 



(1.18) 

+ So H (*) °v (*) + Zuqp (t) tup (t) a pq (t))dt 
This demonstrates the linear microduality principle in the following specified form 
of the table (pT23) 



Filtering q 


A — /Li -1 


a 


e 


J T K 


GJ 


SJ 


Control p 


A - f i- L 


a 


7 


L T 


JH 


m 



showing the complete symmetry under the time reversal and exchange of (q,p), 
in which the coordinate observation is seen as completely dual to the feedback of 
momentum. 



2. Quantum Dynamics with Trajectories 

This section highlights the differences between quantum and classical systems 
and introduces the problem of quantum observation and its solution in the frame- 
work of open dynamics. In orthodox quantum mechanics which treats only closed 
quantum dynamics without observations, there is no such problem. However, it is 
meaningless to consider quantum feedback control without solution of this problem. 
After the appropriate setting of quantum mechanics with observation is given, the 
measurement problem is then restated as a statistical problem of quantum causality 
which can be resolved by optimal dynamical estimation on the output of an open 
quantum system called quantum filtering. 
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Quantum physics which deals with the unavoidable random nature of the mi- 
croworld requires a new, more general, noncommutative theory of stochastic pro- 
cesses than the classical one based on Kolmogorov's axioms. The appropriate 
quantum probability theory was developed through the 70s and 80s by Accardi, 
Belavkin, Gardiner, Holevo, Hudson and Parthasarthy [45l [T7l [19l l4l H6l l47l |48] 
amongst others. 

The essential difference between classical and quantum systems is that classical 
states, including the mixed states, are defined by probability measures not on prop- 
erties but events. This is because the properties of classical systems are described 
by measurable subsets A C CI forming a Boolean cr-algebra 21 on the space of clas- 
sical pure states, the points lo G CI of a phase space. In principle they all can be 
tested simultaneously and identified with the events represented by the indicator 
functions 1a (w) of A £ 21 on the universal observation space CI. They are building 
blocks for classical random variables described by essentially measurable functions 
with respect to a probability measure P on 21. The algebra of all such complex func- 
tions a : CI — > C with pointwise operations is denoted by A, while L p (0,P) with 
p = 1, 2, oo stands for the subspaces of absolutely integrable, square-integrable and 
essentially bounded functions f,g,b 6 A respectively. Note that the Banach space 
M = L°° (Q,P) is a commutative C*-algebra (see Appendix 1.1) of the algebra A 
with involution *:A9ai— >a*GA defined by the complex conjugation a* = a. 
Moreover, it is W*-algebra since M has the preadjoint space = L 1 (CI, P) such 
that M* = M with respect to the standard pairing 



defining the expectation on M* as E [/] = (1, /). 

In quantum world, unfortunately, there are incompatible properties correspond- 
ing to inconsistent but not orthogonal (i.e. not mutually excluding) questions such 
that, if the infimum P A Q is zero, it does not mean that P _L Q. These questions 
cannot be surely answered simultaneously, i.e. tested with simultaneous events on 
any universal measurable space CI, and they cannot be represented in any Boolean 
algebra. Since the incompatibility is measured by noncommutativity of orthopro- 
jectors P and Q representing these questions as Hermitian idempotents on a Hilbcrt 
space TL of quantum vector-states, the algebra B generated by all quantum proper- 
ties must be noncommutative. The set (£>) of all orthoprojectors P G B, called 
property logic of a noncommutative algebra B, clearly extends any eventum logic 
of commuting orthoprojectors injectively representing the Boolean logic 21 by a 
(T-homomorphism E : 21 — > *P (B) such that = I f° r an y measurable 

(T-partition CI = ^2 Aj . Two normal quantum variables are said to be compatible 
if their orthoprojectors commute, and therefore can be represented classically by 
measurable functions on their joint spectrum space CI, however there is no such CI if 
they do not commute. Since there are many incompatible quantum variables, e.g. 
the position and momentum in quantum mechanics, quantum properties cannot be 
identified with any commuting set E (21) representing a Boolean logic 21. 

2.1. Quantum causality and predictions. Almost simultaneously with Kol- 
mogorov's functional formulation of classical probability theory von Neumann |49j 
gave another, more general operator formulation, aiming to lay down the foundation 
of quantum probability theory. It deals with not only commutative W*-algebras, 



(2.1) 
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called von Neumann algebras when they are represented as algebras of operators 
on a Hilbert space Ti with involution as Hermitian conjugation * and unit as the 
identity operator / on Ti. In order to understand the relation between these two 
formulations it is useful to reformulate Kolmogorov's axioms in terms of von Neu- 
mann (wise versa is impossible in the case of noncommutativity of the operator 
algebra). Any random variable a G M can be represented by the diagonal opera- 
tor a of pointwise multiplication ag = ag in the Hilbert space H = L 2 (f2,P) such 
that the abelian (commutative) operator algebra M = {a : a G M} is maximal in 
the algebra B (H) of all bounded operators on H in the sense that M = M'. Here 
M' = {B e 6(H) : M, B = j with [M, B] = |AB - BA : A G M j stands for the 
bounded commutant of M, which obviously coincides on H with the commutant 

(2.2) i'* = {B: [1a,B] = 0,Ag2(} 

of the Boolean algebra 1% = {1a : A G 21} of all diagonal orthoprojectors (the 
multiplications by 1a) generating M. Note that the commutant B = M! of any 
nonmaximal abelian subalgebra (Ti) is a noncommutative W*-algebra with 

strict inclusion of Ai as the center B R B' of B. Thus the simple algebra B = 
B (Ti) is the commutant of the abelian algebra of scalar multipliers M. = CI which 
is generated by the trivial Boolean algebra 21 — {0, il} represented by improper 
orthoprojectors P% — O, Pq = I. The noncommutative algebra A cannot be 
generated by any Boolean algebra of orthoprojectors as the commuting Hermitian 
idempotents P 2 = P = P* in B (Ti). 

Quantum causality, assuming the existence of not only properties but also ob- 
servable events, requires that all quantum properties related to present and future 
at each time-instant t must be compatible with all passed events. This makes an 
allowance for simultaneous predictability of incompatible properties at least in the 
statistical sense. However the usual quantum mechanics, dealing only with irre- 
ducible representations B = B (Ti) of quantum properties but not with the events, 
is causal only for the trivial eventum algebra of improper orthoprojectors {O, /} on 
Ti.. This is why any nontrivial causality requires an extension of the orthodox frame- 
work of quantum mechanics to quantum stochastics unifying it with the framework 
of classical stochastics in a minimalistic way allowing the distinction between the 
future quantum properties and past classical events. This program was completed 
in [H [24] on the basis of quantum nondemolition (QND) principle [HI [20l [27] as 
an algebraic formulation quantum causality. The past events, corresponding to the 
measurable histories A G 2l t ] up to each t G R+ , should be represented in the com- 
mutant of a noncommutative subalgebra By t C B describing the present and future 
on a universal Hilbert space Ti. Thus, instead of a single noncommutative algebra 
B extending the eventum W*-algebra M. generated by E (21) one should consider 
a decreasing family (B[ t ) of reduced subalgebras B[ t C B[ s Vs < t in the corn- 
mutants E (2l t] )' = {B G B : [E (2l t ]) ,B] = 0} of the past eventum logics E (2l t] ) 
representing the consistent histories of increasing probability spaces (f2 4 ],2t t ],P t i) 
in nonmaximal abelian W*-algebras M, t \ generated by E (2l t j). 

The nondemolition principle makes quantum causality irreversible by allowing 
future observations represented by decreasing eventum algebras E (pit) C Bit to 
be incompatible with some present-plus- future questions Q G Bu. Although any 
projectively increasing family of classical probability spaces can be obtained by 
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Kolmogorov construction from a single (0,21, P) with projections k s ] : O — > fi s ] 
inverting the injections K^j 1 (2l s j) C k^ 1 (2l t j) C 21 for all s < t such that P s ] = 
P t j o k^ 1 = Po k^j 1 , however this projective limit may not be compatible with any 
noncommutative algebra By t . Thus the maximal W*-algebra B = B[ , satisfying 
the compatibility condition B' = Ai ] with the initial central algebra A4 ] = BnB', 
coincides with the decomposable algebra Ad'^ which is not compatible with the total 
eventum algebra Ai = VAi t ] except the case Ai = Ado] of absence of innovation 
Ai r ] = Ai t ] for all r and t. The latter with Ai = CI is a standard assumption in 
the orthodox quantum mechanics dealing in the absence of observations with the 
constant By t equal to B(TL). We may assume that Ai ] = CI corresponding to 
trivial initial history, 2l ] = {0,0} with P ] — 1 on a single-point fi ] = {0}j which 
allows £>[ = B (H[o) ■ 

Note that since all operators B G E (2tti) commute with Ait], they are jointly 
decomposable, given in the diagonal representation of Ai t ] by (2l t j, P t j) -essentially 
bounded functions B t : u) i— > B (w) on fl t ] with operator values B (ui) e B (u) on the 
Hilbert components TL (u) of the orthogonal decomposition f® H (u>) P t i (dui) ~ TL 
corresponding to the joint spectral representations 

(2.3) E(A)~ l A (Lu)I(Lu)F t] (dcj)=I t] (A) 

Jn t] 

of commuting orthoprojectors E (A), A G 2l t j. 

Quantum state (See Appendix 1.2) consistent with the trajectory probability 
space (fi,P, 21) is given as the linear positive functional (vj\Q) — (w,Q) by a 
Hermitian-positive w — w* mass-one — 1 operator w h (g B in usual 

or a generalized sense as affiliated to B) defining the probability measure P as the 
projective limit of 

(2.4) P t] (A) = H£(A))=P(A), A eOt] 

(where (w,B) = tr [Bw] for the semifinite B). It is called the (probability) den- 
sity operator for B since it defines the probability Pr [Q] — (ru\Q) G [0, 1] of any 
quantum property described by an orthoprojector Q G A. Since Q G ^ (^[t) is 
compatible with each eventum projector E (A) for A G 2l t j, the property Q is 
statistically predictable with respect to all past events due to the existence of a 
posteriori conditional probability 

(2.5) Pr[Q|A] = — ^HQ25(A)> VA:P(A)^0 

such that Pr [Q] = P (A) Pr [Q\ A] + P (A- 1 ) Pr [Q| A x ] . Note that (zu\QE) is not 
positive and even not real without the compatibility of Q and E. This leads to the 
existence of the posterior quantum states on B[ t given by the conditional expecta- 
tions 

(2.6) e t [B\u} = (w [t (uj)\B(w)) VBeBnM' t] . 

The posterior states are defined as classical stochastic adapted processes u \— > 
U7r t (oj) with density operator values affiliated to the components B(u>) of B\ t , and 
the corresponding conditional expectations defined as positive normal projections 
B n A4jj — > will be denoted as e* = o7* t . 
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Theorem 1. Let w be a normal state on B. Then the induced state vou on the 
relative commutant B] t C BflMJj of the eventum algebra Ait] * s given as classical 

expectation (u7[ t ,B) = En t] w* t (B) in terms of pairing 

(2.7) (m [t \B):=f (w [t (w) \B (w)> P (du>) 

Jn t] 

on Bt = B H A^jj wrai/i the posterior density operators w[t h 23* uniguely defined by 
the conditional expectation e : Bt —> Ai t ] as positive integrable function for almost 
all lj G f2 t i 

Proof. Since w t is normal state on B t , equivalent to the space of essentially bounded 
functions on (O t i , 2l t ] , P t j ) with operator values in S (w) C B(H(lj)), it is uniquely 
defined by an essentially integrable function w t (u>) in (|2.7[) with operator values 
in £>* (cj). Obviously it is the density operator for the posterior states as the con- 
ditional expectations defined on Bt with respect to the central Abelian subalgebra 
Ait] ~ L°° (yi t ],F t ]j by the Radon-Nikodym derivatives 

(2.8) e*[5M:= Ii» BeB t 
K 1 1 ' J a il3 A\H P(A) 

where the limit is understood for and almost all u> € Q t ] m the same way as in the 
classical case. □ 



Note that in the most important "white noise" cases considered in next sec- 
tions, all TL {lj) with oj e Q t ] are isomorphic to a single Hilbert space Tlu of a 
decreasing family \H\t) such that B t (u>) S B \H[t) , uniquely defined for almost all 
w G f2 t j , represents B £ Bt as a bounded operator on TL = H t ] ® TLu commuting 
with I t ] (A) = 1a ® I[t, where H t ] = L 2 (Q t ],F t ]). Representing H as H t ] ® W[ 4 in 
the case 7i (cj) = H[ t G O t i , the posterior states, described by the positive mass 
one operators W] t (w) = w^. in H[t, can be considered as the conditional states on 
the operator algebras B[ t — B \ H[t), controlled by the history trajectory w. Thus 
the above quantum causality setting gives immediately the posterior states Wy t as 
the states for quantum present and future conditioned by the classical past without 
any reference to the projection or other phenomenological reduction postulate of 
quantum measurement. This is main advantage of the extended, event enhanced 
quantum mechanics, or eventum mechanics, which allows treatment of the observ- 
able events on equal basis with other quantum properties of the system. It can 
be shown, see (|7.5p , that any reduction postulate of the operational quantum me- 
chanics can all be derived from QND causality, and this principle is also applicable 
to the continuous measurements in both in time and spectrum where projection 
postulates fails. 

Thus posterior states provide the optimal in mean quadratic sense Bayesian 
estimators for any number of unobservable quantum noncommuting variables x G 
B\H[t) or properties from future, given the observable history u G fi t i. For this 
reason, one can consider quantum measurements in this " nondemolition" setup as 
a form of quantum filtering. 

We now describe an appropriate Markovian model for the time-continuous in- 
teractions between the open quantum system and the field. 
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2.2. Quantum open dynamics and input-output. Quantum Markovian dy- 
namics with observable trajectories, which entered into physics in the 90's in terms 
of stochastic transfer-operators or stochastic Master equations, define the phe- 
nomenological "instruments" of observation without giving any microscopic dy- 
namical model in terms of the fundamental Hamiltonian interactions. In fact such 
approach is equivalent to the earlier operational approach based on the instru- 
mental transfer-measures (See Appendix 1.3), and its starting point corresponds to 
already filtered Markov dynamics in the classical case. Here we describe the general 
scheme for underlying Hamiltonian interaction models with continuous observation 
for open quantum dynamical objects in terms of quantum stochastic evolutions in 
parallel to the classical stochastic models with partial observation, following the 
original Belavkin approach suggested in [4, 5, 24J . 

Let us fix a quantum probability space (Ti., B, w) and an increasing family B s \ C 
B t ] , Vs < t of W*-subalgebras B t \ C B containing the compatible histories E (2l t j ) C 
2?j t and nontrivial present B[ t ] = B t \ H By t for each t, assuming that each B t \ com- 
mutes with future At C £>[ t , B t ] C A' tl being generated by only nonanticipating 
questions Q G A' t fl B. The future is quantum noise which is described by W*- 
subalgebras At Q B, forming a decreasing family At C At+ S Vi, s > with trivial 
intersection such that we may assume that VSj] = B. Moreover, we shall assume 
that the family {B t \) as well as (M t \) form the W*-product systems in the sense of 
W*-isomorphisms 

(2.9) B t] ®At ~ B t+s] , M t] ®M s t ~ M t+S] , 

where A\ = At H B t+S ], Mf = M t l~l M t + S ] and (M t ) is decreasing family of W*- 
algebras M ( C B[ t generated by future events E (2lt). This implies that the family 
(Af) satisfies the product condition such that Ao ~ A t ^)At^)At+ s for any t and 
s > 0, and similar for (.Mf), corresponding to the split property 

(2.10) n = n* x Q r t x n s t+r , 21 = % ® % ® %t+ r 

of the measurable trajectory space. 

Quantum open object under the observation is represented at each time t by a 
past-future boundary W*-subalgebra bt C B[ t ] such that it is adapted with respect 
to the family (£>*]) quantum stochastic process (in the genral sense [3]), nonantici- 
pating futures (At) and satisfying causality condition with respect to the histories 
(A4 t ])- We may assume that each bt represents a fixed b or a variable boundary 
W*-algebra b (t) by a W*-homomorphism Tr t of b (t) onto bt, with b (t) taken in the 
initial algebra 23[ ] , say. Due to the causality condition the product 

lit (A, q) = E (A) 7T t (q) VA G 2l 4] 

defines for each t an adapted transfer-measure n t (A) : b (t) — > B t ] (see Appendix 
1.4) with W*-homomorphic values normalized to the history eventum projectors 
E(A). Obviously W*-algebras b t and Al are in Bt = B t+S \ n B[ t , as well as b t+s 
and W*-algebras M s t . 

Following [4] we shall say that quantum open object b (t) with eventum history 
E (2ltj) is dynamical with respect to (At) if 

M s t V b t+s C b t V A s t Vi, s > 0. 

This is equivalent [4] to the existence of quantum flow with observations described 
as follows on the co-images b (t) = 7r' (b t ) of the boundary algebras b t . Assuming 
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that the W*-algebras Bt are generated by bt and A s t: we can always consider the 
dynamical quantum open object with bt = Buy 

Theorem 2. Let tt : bt — > b (t) be normal injections inverted by the dynamical 
representations n t , and let B t yM.t\ form the product systems i2.9\) . Then there 
exists a transitional spectral measure 

(2.H) T t r (A,q)^E t (A)j t r (q) 7 qeb(t + r) 

on 2l[ with values in B? t — b (t) >S>Al given by adapted a-homomorphisms E : 2l t — > 
b (t) ®At and a Heisenberg flow (7* ) of causal tensor- adapted W*-homomorphisms 
7* : b (t + r) -> E 1 (%)' n B\ t such that 

(2.12) 7*°7' s +r = lUs Vr, S >0, 

(2.13) 7 *- r (S t (A)) = E^ r {A) Vt>r 

under the trivial extensions onto Bu = b (t) ®At ■ 

Proof. The representations n t as well as 7r* _r can be trivially extended to the 
adapted W*-homomorphisms with respect to the identity maps id respectively on 
At and At- r by virtue of commutativity B t ] C A' t as n t (q A t ) = lit (<?) A t and 
Tt t ~ r (qt- r A t -r) = n t ~ r (qt-r) ® A t -r respectively for all q £ b (t) and q t £ b t , 
At £ At- This defines the compositions r f~ r = 7r*~ r o Tt t of thus extended W*- 
representations as tensor- adapted W*-homomorphisms r ) t r ~ r : Bu — > Bu—r on 23[ t = 
b (i) (8) -At trivially extending the b (t) — > b (i) ® -4[_ r and satisfying the hemigroup 
condition (|2 . 1 2[) such that 7q = id (J3u) for each t. Obviously these extensions 
satisfy causality condition 

n t (B [t ) c (a t] )' , 7 * (B [t +r) c (an' , 

where the eventum projectors 23* (A) £ are defined for any A £ 2l t as 7r* (_E (A)) 
by the extended injections 7r* : b t V A — > b (t) ® «4t as right inverse of the extended 
TT t on E — TT t (E t ) £ E(% t )- The second condition (|2.13[) simply follows from 
7r f o 7* = 7r t+r due to 7r t _ r (£'*) = £ for any r £ [0,t] and E £ £(2tt). Thus 
the QS flow with nondemolition observations can be described in terms of the 
homomorphic transitional measures (|2.11[) with (|2.12p and (|2 . 13[) satisfying the 
hemigroup composition low 

(2-14) (A[_ r , Tl (A s t ,q)) = TlZ (A[± r s , q) 

where A r t +^ = A[_ r xA[e 2l^ r s and q £ b (t + s). □ 

Corollary 1. The dynamical QS object is Markovian in the usual sense |3j if the 
initial state vj = n7[ on B — £>[ is product state w ~ w t \ ® g t for any t such that 

(w\B t] A s t ) = (w\B t] ) (g t \Al) . 

It is operationally described is such state by the hemigroup of reduced transitional 
measures 

(2.15) 7?(A,q) = g* t [? t r (A,q)], 
where : B[q — > B t \ is conditional expectation defined as 

(2.16) <tJ7[ |S t] ) = (w t] \g* [B t] ]) Vw t]) B t] £ B t] . 
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They satisfy the operational Chapman- Kolmogorov equation 
(2-17) T?~ r (A[_ r , Tl (A!, q)) = T r *+J (AJ+', q) 

as a normal completely positive map b (t + s) — > b (t — r) for each product A[+* of 
A r t _ r €W t _ r and At e Of. 

Remark 1. TTie event representations E = 7r* (E) are usually given by input 
a-homomorphisms I : 2l| — > .4£, ^(A) = i (1a) as E* (A) — v f (I (A)) m ierms 
o/ a two side adapted W* -representation i : Mf — > .4* and a hemigroup (u*) o/ 
interaction isomorphisms t>* : fa (i + r) — > b (t) <§>.4£ s«c/i £/iat u'|fa(t + r) = 
7* . ffere w* defines the output representation induced on an input eventum algebra 
t(M t ) C A /or M t = L°° (fit, Pi) 6j/ the limit v t = ]hn r v t r \A t of v\\A\ which is 
well-defined on each At = \lrA\ due to the localization property t>*_|_ r |*4£ = v r \A r t 
for allr,s > 0. Note that the localization property simply follows from the hemigroup 
condition and the normalization v l (Io) = h for these v\, extended adaptively also 
on A G Aq and A t + r € At+ r such that 

vl (Al <g> B ® At+r) = A f <g> vl (B) ® At+r- 

The quantum free evolution is usually described by a semigroup (6 S ) of endomor- 
phisms 6 S : £>[ — > B[ shifting isomorphically any A t onto A t + S with trivial action 
on fa. QS Hcisenbcrg flow (7*) with observation over a constant algebra fa (t) = fa 
is called covariant with respect to a shift semigroup (0 S ) acting also on 2l by the 
shift of any 2l t onto 2l t+s = 6 S (21*) if 

(2.18) 7^°0t+ s = O t o0 8 , E t o6 t = 9 t oE , 

where ^ s = 7° o S and 6 t is extended on the W*-algebra B\q = fa ® Aq by 
0t (9 ® -<4o) = 9 ® #t (-4.0 )■ This defines a Heisenberg 6>-cocycle 7 S = 7° corre- 
sponding to the semigroup (i? s ) of W*-endomorphisms ^ s = 7 S o S of the algebra 
B, satisfying causality condition d s (B) C E 1 (2l s j) . Note that the shift semigroup 
can be extended to a group {8 t : t € M} on = .4o<S>b<S>.Ao, where *4o is an inde- 
pendent copy of the algebra Ao, with 9t transforming each segment AJ. onto A T S 
for any positive T,r,s and t = r + r + s, shifting At onto .4o similar to the in- 
verted shift 6-t ■ At — > Ao, with backward transformation of each Al onto A T r for 
r+r+s = — t. This free group dynamics 8 t defines the reversible quantum dynamics 
interaction on such B by one-parametric group of i? t = v t o t+s extending $ t onto 
£> by interaction W*-automorphisms u t = u° on B[ = b®Ao, acting identically on 
Ao for any t > 0, with the identical action on _4o and the reflected cocycle action 
v t = _t o v-t o 8 t + s on B ] = Ao ® fa for t < 0, where u s is defined on B ] exactly as 
v^ 1 on B[ for any s > 0. However the reversible quantum dynamics on such non- 
commutative B cannot satisfy the causality in both directions of time with respect 
to a nontrivial eventum algebra E (21), except the case of absence innovation as in 
the conservative quantum mechanics without observation. To keep the causality in 
the positive direction of time one must replace the nonabelian .4,0 by the smaller, 
abelian subalgebra M , a copy of the eventum algebra M = / (Mo), which makes 
6 S and ^ s irreversible on B = Mo®b®Ao with noncommutative Ao- 

3. Quantum Stochastics of Eventum Mechanics 

In this section we consider quantum noise models defining quantum Markov- 
ian dynamics with continuous nondemolition observation and show the quantum 
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filtering equations derived from this models. Such observation can be based only 
on indirect measurement of a quantum open object via a coupled channel repre- 
senting a classical measured process m in a bath Aq which is usually assumed to 
be initially independent of the quantum object b. We shall consider the measure- 
ment processes having initially independent increments mf = rn +s — = m (I s ) 
as random measure on the intervals I s = [t, t + s). They generate independent 
W*-algebras Mf = L§f(f2£,Pf) = M(If), and an input quantum process with 
the increments M t s = t(raf) = M (If) generates independent eventum algebras 
i(M|) C B (T£) on a Hilbert space To with respect to an initial unit state vector 
Xo satisfying the divisibility condition 

(3-1) T ~ Tl®T 8 t ®T t +s, Xo~Xo®xt® Xt+s, 

such that it induces the initial product probability measure Po = P ® Pf ® P*+ s on 
the measurable space of observable trajectories under the split condition (|2 . 10[) . 

An appropriate candidate for such Hilbert space suitable to accommodate any 
kind of classical independent increment process is Guichardet-Fock space To = 
r (So) over the Hilbert space So of L 2 -functions £ : t i— > t on = {t > 0} with 
values in a Hilbert space t such that Tl — T (£ t s ) with £ t s = L 2 (I s — > 6) (see the 
definitions in [24| . [50] summarized in the Appendix2). There are sufficiently many 
product vectors in J-q, called exponential vectors {£® : £ S So}, generating Fock 
space To such that any coherent state vector 

(3-2) X e = e-««£® Uf £ = / U(t)\\ldt 

Jo 

defines a product state (v^l^) = (x^l^X^) on an y subalgebra .4o Q B(To) satis- 
fying the divisibility condition Ao ~ A* <£> A s ® At+ S for any t, s > such that 
each A s = A (I s ) is represented in B (T s ). However there is only one shift-invariant 
such state which is given by the vacuum vector \o = ^0 corresponding to £ = 0. In 
fact, not only shift invariant but any infinitely divisible normal state on Aq can be 
induced by the vacuum state ip on B (To) by choosing in general time dependent 
Hilbert space t (t) in a canonical way [51], [52]. In particular, since the state on 
the abelian algebra Mo defined by the probability Po of any classical process with 
independent increments is infinitely divisible, any such state can be induced from 
the quantum vacuum state by restricting it to the abelian part t (Mo) C Ao given 
by an adapted input W*-representation l : Mo — * .Ao such that i (Mf ) c A s t for any 
t,s > 0. It defines the infinitely divisible probability measure Po as 

(3.3) P t s (A):=(^|/(A)^)^(^|/(A)) VA e 21*, 

where / (A) = i (1a), by the vacuum vectors Xt — <^0 01 ■ 

There are two basic processes with additive independent increments which can be 
realized in Fock space with finite-dimensional i — C d : Wiener vector- valued process 
= (wl) ieI r which we index by a subset Iw Cj {1, • ■ ■ , d}, and Poisson compound 
process n (I ) = (n') igJ which we index by another subset In Q {l,...,d} (It 
should be thought as diagonal matrix- valued rather than vector). Their differential 
increments dm = m (dt) satisfy quite different Ito multiplication tables 



(3.4) 



dwidwk = Sikdt, dnidn k = 8jS J k drij, 
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where the summation rule is applied over j. Their canonical input representations 
M* = i (m\) in Fock space over £q = C d ® L 2 (iR^) are defined as 

(3.5) Wl = Af(t)+At(t) = $t[A+(t)], 

(3.6) N* = A*(t)+A+{t)+AL{t)+A±{t) 

in terms of four basic operator- processes in J-q defined in Appendix B. These are cre- 
ation A+ (t) := A+ (l£) (row- valued, A+ = (Af)), annihilation A°_ (t) := A°_ (l£) 
(column- valued, A°_ := (A k _)), exchange A°(t) := A° (I§) (matrix-valued, A° = 
(Af) ) and preservation A+ (t) := A+ (Iq) (scalar- valued, A~t (t) = tl). These 
canonical processes forming a pseudo-Hermitian matrix A = (^4^)^T°'^ = A* under 
the involution (At L )* = (A L * K ) with respect to the index reflection — (— ,°, +) = 
(+, o, — ), satisfy pseudo-Poisson multiplication table 

(3.7) d^d^ = ^d^ \/L,ve{i,+},ii,K€{-,k} 

of quantum stochastic calculus discovered as a noncommutative generalization of 
the classical Ito-Poisson table by Belavkin in 24J. Note that from (|3.7p it follows 
that 

dWidNk = SldAl + 6 lk dt, dN,dW k = S l k dA+ + S tk dt 

which cannot be realized in the classical category of commutative processes in which 
it always dwidn k = = dn k dwi. Thus the joint operator representation of two 
types of basic classical processes with independent increments is possible only in 
splitted Fock spaces such that Iw H Jjv = corresponding to orthogonal subspaces 
iw -L %n, which reflects the classical split 

n = n w x n N , 21 = a w ® % N , p = p w ® Ptv- 

The four basic processes A v ^ form a linear basis of quantum Ito ★-algebra as 
noncommutative integrators for the increments 

(3.8) Ml = J2 K " M dA » = l * ^ 

defined by four integrants K — (K^)^~~ , ° as operator- valued functions integrable 
in a quantum-stochastic [24] . These quantum stochastic integrals satisfy ^-property 
such that if (K)* = if (K*) under the involution (Kt v )* = (K^), and the Ito 
product rule 

(3.9) d(M*M) = d (M*) M + M* (dM) + (dM*) d(M), 

where dM = K^dA 1 ^ (the usual summation convention is assumed), with the Ito 
correction calculated as 

(3.10) d%(K)*d%(K) = d%(K*K) 

for adapted quantum stochastic integrators K£ (t). 

Note that in the case of a single degree of freedom d — 1 this quantum Ito table 
reads in the Hudson-Parthasarathy (HP) form [IB] as 

dA°_dAt = dtl, dA°_dA°=dA°_, 
dA°dA° = dA°, dA°dAt =dA+, 
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with all other increment multiplications vanishing, in terms of four scalar-operator 
processes A£ with pseudo-Hermitian property written as 

At (t) = At (t)* , A+ (t) = A°_ (t)* , A° (i) = A° (<)* . 

3.1. Quantum mechanics with observations. Interaction automorphic evolu- 
tions (v r ) r>0 on the the tensor product £>[ — b ® -4o 01 simple algebras b = B (h) 
and Ao = B (T ) are usually described by right unitary cocycles {U r : r > 0} as 
v t (B) — U t BUt . The cocycles satisfy the operator hemigroup condition U*Ul +r — 
U* +s similar to (|2.12[) in terms of the shifted operators U* = 6 S (U r ), where (0 S ) is 
the semigroup of right shift W*-endomorphisms 9 S : Ao — > A s C „4o describing free 
evolution of the bath by trivial placing of each subalgebra A (l?_ s ) onto A (If) as 
A% = B Hudson and Parthasarathy [48], [50] derived a QS forward equation 

(3.11) dU t = U t (RdA° + R+dAt + R~dA°_ + R+dt) 

for the unitary cocycles in Fock space To = T (L 2 (Rj)) defined as the solution 
of QS integral equation Ut = Uq + i (UR) with four adapted operator-valued 
coefficients R^ and Uq = I. They gave the necessary conditions for unitarity of this 
solution written in terms S = I + R as 

(3.12) S* = S-\ R- = -R* + S, 29? (R+) = -~R* + R+, 

where $t(A) denotes the Hermitian part (A + A*) /2 of an operator A (the suffi- 
ciency was shown only for the constant bounded initial- valued coefficients R^)- 

It is important for control to have also the sufficient unitarity conditions in 
the case of time-dependent and multi-dimensional noise A. As it was proved in [55] 
under the natural QS-integrability conditions, the relations f|3. 12[) are also sufficient 
for the uniqueness of unitary solution even in the case of time-dependent infinite 
dimensional coefficients with values R$ (t) = (s(j (t) — 5^1)®/ in the initial operator 
algebra b. In terms of these s£J (t) the multidimensional version of HP equation can 
be simply written in Belavkin's ★-algebraic notations [Ml [SS] as 

(3.13) dU t = U t (t) - 5tl) Al (dt) 

where usual summation convention over all /i, v G {— , °, +} can be restricted to 
the domain [i < v under the order — < o < + of the the triangular matrix 
s = (s£)^~Z'o't with zero operator entries for /i > v and §Z = 1 = §t (the 
usual summation convention then can be applied). In this notations the algebraic 
relations between nonzero matrix elements s£J generalizing the unitarity conditions 
(|3.12p are simply expressed as the pseudo- unitarity s* = s^ 1 of the operator matrix 
s in terms of the pseudo-Hermitian adjoint matrix s* = (sZ„*j. 

From the quantum ltd rule (|3 . 9[) applied to Heisenberg QS flow X (t) = j t (x) 
given by interaction dynamics as j t (x) = Ut(x®I)Ut = vt(x®I), and the quantum 
Ito multiplication table (|3.7p , we obtain the general QS Langevin equation 

(3.14) dX = (EJJ (X) - X5i) dA^ = (E (X) - XI) ■ dA, 

Here E(t,A) = S (t) (X ® 1) S (t)*, called QS germ, is given on X = X (t) by 
matrix-function E (t) which is defined as a triangular matrix (E^)^~_'°'^ of six 
time evolved maps 

E£ (t, U t (4 ® I) Ut) = U t (a£ (t, x) ® /) U* t 
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given on b by er (x) = s (i 1) I*, where 1 = (££)£ZZ'o'+) with trivial aZ = 
id (b) = a\ and (x) = for /i > v such that the summation convention can be 
applied only for \x < v. It was proved in [53j under the natural QS-integrability 
conditions that the unitality and *-multiplicativity of the germ 

(3.15) <r (t, x*x) — <x (t, x) it (t, x) , cr (t, 1J = 1, 

are the necessary and sufficient conditions for the existence and uniqueness of the 
unital *-homomorphic solutions X (r + s) — 7^ (x) to the Langevin equation (|3. 14|) 
with X(t) = x. 

The composition f £_ r (x) = g* [Y t _ r (x)l with noise conditional expectation 
(|2.16|) . defined by the vacuum state g r — ip , describes a dynamical hemigroup 
(fg) (or semigroup (f s ) in the stationary case) of unital completely positive maps 
f r s : b — > b on operator algebra b C £>(()). This bath expectation, given by so- 
lutions X (r + s) = j 7 s (x) of QS flow equation ()3.14j) with X (r) = x <£> I for the 
evolved on time interval I* operators x G b, satisfies the master equation [561 157] 
-riTt (x) — f t (A (x)J with a Lindblad generator A which can be written as a linear 
conditionally positive map in x £ b in the form 

(3.16) X(x) =J2 K " & k( ± ) Kk + L± + ±L *- 

Here K l — —R l + , L = RZ satisfies the condition L + L* = — s 22 li K l *K l and 
a\ (x) = Y^j S'jxSj* (R\. — §Z, RZ = sZ, Sj = to denote that these operators 
belong or affiliated to the algebra b). In fact, for quantum coherent control we need 
a time dependent version of this equation in the following decomposed form. 

Theorem 3. Let (p*^ be coherent conditional expectation B[ — ► £> r ] defined on 
B[ = b <£> Aq as normal positive projection such that 

{*>r]\>P%,t [B]) = (V r] ®X( r \B\il> r] ®X( r ) 

for any vector product-state vj t ^ = ® (pgr with r\ G f) and coherent vector AS.ty) 
in Tq. Then the dynamical map f (t, x) — ip* ^ [Tt_ r (&)] = r£_ r (£; x) satisfies 
forward evolution equation 

(3.17) ±f(t,i)=f {t,X(t,x)), f(r)=Id(b) 
with Lindblad type generator decomposed as 

(3.18) A (i) = i [H, x] + \ R (x) + K°* (<r° (x) - xS°) K° 

where H (t) = TiQ [(2£ (t)* R° + + R~)] , K° (t) =£(t)- R° + , a° (x) = S°xS°* and 

i 

Proof. Indeed, it can be shown [58j that conditional coherent expectations f T t _ r (£, x), 
evaluated from the QS flow equation (|3.14j) . satisfy Lindblad type equation with 
generator 

a (t) = az + tot (t) + e (*)* *° + + e (ty (*) 

where p\ (x) = a\ (x) — x8\. Using HP conditions (|3.12p in multidimensional form 

QO* qO j~0-. nO* Q° D— D * DO p — * p — 

o o — o 5 I o — Jt-Q . -'*'-)- — I Jl 1 
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as pseudo-unitarity conditions T_ v — in terms of operator matrix elements of 
inverse T = S^ 1 to the triangular operator matrix S for at u (x) — YliS?%S-*i 
this generator can be written as 

\[x}= J2 (sz,+esi L )x(sr+es:)*-e^ 

i=-,o,+ 

where £ (t)* £ (t) = ||£ This gives Q3.18P after taking into account again the 

unitarity conditions. □ 

Given a normal quantum state <j (r) on b at time r, the hemigroup (f£) defines an 
averaged coherent controlled non-Hamiltonian state evolution (t) : t > r) of the 
quantum dynamical object by composing it with the adjoint CP maps r* = fj*_ r 
as q (t) — t\ (f, (t — r)). It satisfies the hemigroup master equation 

(3.19) ^s + Ks + sK* =^L\U*, 

i 

where the operators K = — L* and L 1 — are defined by left adjoints L* = 
L\L^ = L\ to 

l = r- + + CK - \d, U = r; + es° 

with respect to the (b*|b)-pairing: (Lh\x) = (c;\Lx) (which are usual Hilbert space 
adjoints, L+ = L* , Lj* = L* in the case of the trace pairing (<;\x) = tr [?*x]). This 
master equation is usually written in the Lindbladian form as 4r<; = A (t, q), and it 
is a particular case of the general QS Master equation (|7.13|) derived in [59"].[6"D]. 

3.2. Quantum nonlinear filtering equations. A time continuous measurement 
of each Wiener process w\ represented by the field quadratures (|3.5p after interac- 
tion with quantum object as Yf = UfWfU^, due to locality for any t' > t, realizes 
an indirect measurement of the evolved generalized coordinate Qi (t) — 29? [Li (t)}, 
where U (t) = U t (§i <Z> I) U?. This can be seen from the quantum ltd formula 
(13~§j) , (j3~TU|) applied to the output operators v t {Wj) = UtWfUf: 

(3.20) dv t (Wl) = 23? (L* (t) dt + dA+) = Q, (t) dt + dW* 

Similarly the output process corresponding to the field counting process (|3.6[) as 
Yf = Ut'NfU*, for any t' > t is given by v t {N\) = U t NlU% as the QS integral of 

(3.21) dv t {Nf) = Li (t) L* (t) dt + 25ft (L* (t) dA+) + dA\. 

Classically, filtering equations are used when we need to estimate the value of 
dynamical variables about which we have incomplete knowledge due to an indi- 
rect observation. For example, the Kalman-Bucy filter [61], [62] gives a continuous 
least-squares estimator for a Gaussian classical random variable with linear dy- 
namics when we only have access to a correlated, noisy output signal. Since closed 
quantum systems are fundamentally unobservable (hidden) unless they are open, 
e.g. disturbed by quantum noise processes (c.f. (|3. 14|) such that equations (13.201) 1 
and (|3.2ip have nontrivial input from the quantum object in terms of the non- 
Hamiltonian part L l = s* of the Lindblad generator, filtering of quantum noise 
plays an important role in quantum measurement. As it follows immediately from 
the localization property of quantum interaction evolution due to the hemigroup 
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property the output operators are self non-demolition (i.e. mutually compatible 
at all times) and satisfy the quantum non-demolition (QND) condition 

(3.22) [X(s),Yf]=0 Vs < t,i G 

with respect to any evolved quantum object process X (s) = v s (x <g> 7) = 7 S (x). 
Belavkin was the first to realize that an optimal estimation without further dis- 
turbance is possible in the general quantum open dynamical models when based 
on any output QND measurements [18] , [44] , [20] , [5] . He constructed the quantum 
filtering equation which describes the evolution of the optimal estimate given by 
the density matrix conditioned on a classical output of the noisy quantum channel. 
This is used to estimate arbitrary object variable X (t) h bt which are driven by 
environmental quantum noises. The QND condition insists that the expectation of 
X (i) is not disturbed when we measure Yf for s < t. As it was already pointed 
out by the Section 1 Theorem, this is necessary and sufficient for the existence of 
a well defined conditional expectation of X (t) with respect to past measurement 
results of Y^ . 

Let M t \ be the the history abelian W*-algebra W t ] C B (Ft]) generated by the 
output operators {v r (W[) : r G I , i G Iw} for an index subset Iw Q {1, • • • , d}, or 
another abelian W*-algebra Af t \ C B (Ft]) generated by [v r (Nf) : r G I ,i G In} 
for the same or another index subset In, or the product algebra A4t] ~ Wt]<S)A/j] 
generated by Y^} — {Yf : t G Iq, i G I e } corresponding to the union I e C {1, . . . , d} 
of disjoint index subsets Iw and In- Also let 2?[ t denote the future nonabelian 
W*-algebra generated by the system operators X (s) h 7 S (b) for s > t. From 
the compatibility of the output operators, we have quantum causality condition 
M.t] C $[ t , so again we have a unique well defined conditional expectation e = nj* t 
given by the posterior states rorj on Bn onto A^tj. 

The conditional expectation x* = e*[X(t + r)] gives for any r > the best 
prediction of X (t + r) = v\ (X (t)) as least squares estimator of any operator 
X (t) = Vt (x ® I) evolved to the time t' = t + r conditional on the output operators 
Ye ] and so is equivalent to a classical random variable on the space of measurement 
trajectories f2 t ] — {cj[ : r G I ,i G / e } such that w[ is an eigenvalue of Yf. This 
conditional expectation for is most conveniently written in the Schrodinger picture 
as £' = (^'|t* (£)) in terms of the expected CP hemigroup r* and posterior states 

is defined by the relation 

(?\x) := e*[w t (x ® J)] = (n7[ t |fi (x ® /)) . 

In the case of product state w\ r = <^ ® Q r with infinitely divisible g r realized on 
Ar = B (F r ) by the vacuum state, the posterior state is given for any t > r as 
£ (i) = 0^,~ r (<r) by a hemigroup of nonlinear transformations of a starting 

state <r = ? (r) on the object algebra b, resolving the quantum nonlinear filtering 
(Belavkin) equation 

(3.23) d? (t) = A [?] (t) di + ^ 5' 1 (?) (t) Af, (dt) . 

Here A is Lindblad generator of the adjoint equation (|3.19p . and quantum filtering 
coefficients S l against the innovation martingales 

Mi Q-r) = Yi (I*) — /"" *e*[y<(di)] 
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on the time intervals If where first specified in [63] as the functionals of the posterior 
quantum states ?* = ?,■ (r + s) — <f)^. (?) resolving Ito equation (|3.23|) for ? (r) = ?. 

We present here two separate cases of Belavkin quantum filtering equation cor- 
responding to the diffusive and counting measurements. For rigorous derivation see 
[36] . and for the most general mixed case we refer to [24 . 

3.2.1. Diffusive Belavkin equation. The diffusive version L 63J , 32J, 36J of Belavkin 
quantum filtering equation corresponding to continuous observation of the diffusive 
row-vector Y l w = v t (W^) = (Y?) ieI indexed by a fixed set I e = I w specifying 
the estimation channels is a classical non-linear stochastic differential equation 
given by 

(3.24) d? = A [?] dt+J2 PtfXXi ( dt ) - + L*) dt). 

Here A (?) = £V L\L"-Kq-qK* is defined by the left adjoints L 1 = h\ , K = -L* 
of the Lindblad operators in (|3.16[) . and 

5 l (?) = ?L** + L\ - (? \Li + L|) ? 

is the nonlinear fluctuation coefficient such that (<T (?) |l) = for any ? £ b* with 
respect to the pairing of b and b*. 

3.2.2. Counting Belavkin equation. The counting, or quantum jump version |63[ 
134] . [36j of Belavkin quantum filtering equation corresponding to counting obser- 
vation of the number processes Y l N — v t [Njj — iYl) ieI indexed by I e = In is 
given as a stochastic differential equation in the classical Ito form by 

(3.25) d? = A [?] dt+J2 sl (0 i Y i ( dt ) ~ (Z\LiL*)dt) 

i£l N 

for a counting measurement in the field. Here A (?) = J2i L l qL l * — Kq — qK* such 
that (A (?) \V) = for any ? £ b* and 

(3.26) <f(?) = « l (?)-?, a l (?) L ' a '' 



is the non-linear normalized difference coefficient defined by L % = L\ (=L* for the 
trace pairing of b with b*) such that (5 l (?) |l) = for any ? £ b* with respect to 
the standard pairing of b* and b. 



4. Optimal Quantum Feedback Control 

We now couple the system to a control force (row- vector) u = (ui) ieIf via 
forward (feedback) channels indexed by a finite set If, \If\ = df. The force perturbs 
open quantum dynamics described by unitary cocycle {U r } by making it causally 
dependent on each time interval [r, r + s) on the control segment u s r — u(t) : t £ 
[r,r + s j) such that (r + s) — J7J («*). The family {Z7J (u*)} is assumed to 
satisfy QS equation (|3.13p with controlled pseudo- unitary germ § (t, u) = § (u (t)) 
such that its solution forms a hemigroup 

U t r (u r t )U t + r (u! +r )=U t r+s K +s ). 
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for any composed segment u r t +s — (u[,u r t +s ). This defines a hemigroup (7,) of 
time-dependent interaction dynamics 

7 ; (t,x) = U r u (t)(x ® T)U?(t) - Yt-r «~ r ,i) 

satisfying the controlled Langevin equation (|3. 14|) with S M (t) given by time evolved 
germ & u (t) =& (u{t)). 

Following the original Belavkin's formulation [18j , [3] , [6] of quantum optimal con- 
trol theory, we assume that the quality of a control process on a quantum object 
over a finite period [r, T) with starting product state w\ r — q ® g r on the object 
plus noise algebra b ® A r is judged by the integral expectation 

(4.1) J r (c,«r) = / (w [r |C('u r ,t)>di + (a7 [r |5(w r ,r)) 

of the operator-valued cost junctionals of u r given by the evolved object operator- 
valued measurable positive cost function c (t) : U — > b and a terminal positive cost 
operator s G b in 

(4.2) C (u r , t) = y^l (c(u)) , S (u r , T) = (§) 

for self-adjoint positive system operators c(u(t)),s h A. An alternative problem 
of risk-sensitive control has also been studied by James |11[ 112] where the cost is 
exponentiated to enforce higher penalties for undesirable behavior. 

Coherent control of open quantum dynamics uses field channels indexed by a 
subset If C {1, . . . , d}. It is realized by controlling quantum noise in these chan- 
nels by u e R d f via their coherent states. One can start with uncontrolled dy- 
namics described by Hamiltonian and Lindbladian operators Hq = ^(i?^) and 
Kq = — R\_ for R L + = s L + , l — —,i, defining QS unitary evolution by (|3. with 
arbitrary unitary scattering 5° — s°, and apply coherent conditional expectation 
(p* t to corresponding QS flow y r s (x) with ^ (t) — jrU T (t) defining controlling field 
expectations 

<p U |fiS* [4 (!*)]>= [ S Ui(t + s)ds Vie If. 
Jo 

This effectively results in change of the Hamiltonian Ho and all operators Kq with 
i G // to 

(4.3) H U = H Q + u t Wi (K$ , Kl = Kq + l -u h 

and no change for other Kq with i ^ If (The summation is taken only over i G 
If.) The resulting conditioned dynamics f r u (t, x) satisfies time dependent Lindblad 
equation which can be written in the form (|3.18p as 

(4.4) X u (x) = I [H 2u) x] + X R (x) + J2 K uPl (*) K t 

i,kel f 

with p l k (x) = a l k (x) — S k x, without change of the part \r (but with doubled u in 
the Hamiltonian H u ). 

We are going to consider quantum feedback control problem in which it is nat- 
ural to assume that the forward feedback control channels are disjoint to the set 
I e = Iw U In of estimation. This is achieved by considering coherent controls 
in the channels // such that // n I e = 0. In this case the controlling amplitude 
£ (t) G tf is orthogonal to the subspace t e = tw © tjv of observation channels, so 
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the output equations (|3.20|) . (|3.21|) are not affected by the coherent control which 
will simplify optimal feedback control problem which we solve by applying dynam- 
ical programming to coherent controlled quantum states. The controlled posterior 
density operator £ (t) = can then be obtained from the relevant uncontrolled 
filtering equation by replacing Lindblad generator A = A* in (|3.23p by time depen- 
dent \ u (t) — A (u (t)), and so we now have a controlled time dependent nonlinear 
filtering dynamics u — ^ satisfying Belavkin equation of the form, say (|3.24|) 
or (|3.25p . in which the fluctuating part under the above coherent control assumption 
is independent of u (t) . 

4.1. Dynamical programming of quantum states. In the search for optimal 
control inputs, it is desirable to allow the control to be determined in terms of 
measurement results on the system, particularly in the quantum setting where 
quantum noises introduce an inevitable stochastic nature. A feedback strategy x^o 
consists of measurable maps h (t) which give for each < t < T an operator- valued 
control law u (t) = x{t, Y*J ) as a function of the current and previous commuting 
output operators Y*} = {Y r e : r 6 (0,i]}. Thus the control law x(t) is realized by 
an adapted random vector variable x{t,u> t ^) on the probability space (f2 t ],2l t ],P t j) 
of output measurement results in the value space U of admissible control inputs in 
the spectral representation 

u(t) = f x(t,w t ])Eo (du t] ) = x(t,Y^). 
Jn t] 

We denote the space of admissible operator- valued feedback controls on the interval 
[r,r + s) by Vf,(Y e ). Note that no measurement results are available initially, 
so the initial control it (0) is deterministic and also no controls are applied at 
the termination time T . It is too restrictive to consider only continuous sample 
paths {x(t, a;*) : t > 0}, since for example the Poisson process {3^} certainly does 
not have continuous sample paths. Instead we give the following definition of an 
admissible strategy. 

Definition 1. An admissible feedback control strategy x r = {>c(t) : r < t < T} 
determines randomized control laws u(t) = xlt^Y 1 }) at each time t > which 
realize values in U and form cadlag sample paths {x(t,uj) : t > r}Q. Moreover, an 
admissible strategy x° t shall be called optimal if it realizes the infimum 

(4.5) S(r,?)= inf J r (?,x r ) = J r (c,x°) 

over the space IC r of admissible feedback control strategies, where J r (<;,>< r ) is the 
expected cost for the control process determined by the feedback strategy >c r . 

It is a simple exercise to show that under a feedback strategy, the output oper- 
ators once again form a QND measurement with respect to the controlled system 
operators which justifies the existence of the conditional expectation 

(4-6) e*bC(*]] HO) 

with respect to the output operators Y 1 } for t > r . It is given by the posterior state 
*>r~u = ^r,u (t) a s solution of Belavkin equation (|3.23|) which now has dependence on 

^Right-continuous paths (lim^_,o A t t+h( a; * + ' 1 ) = Mt( a; *)) having well denned left limits 
Mt<V)- = lim h _ Mt-/ l (^ t "' 1 ) 
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the chosen control u r s inputs through the dynamics 7* (u r , x), and, given an initial 
condition q r>u (r) = q, it does not really depend on Y T } due to Markovianity of the 
process ^ u proved in [3] . The existence of this conditional expectation permits the 
following theorem which lies at the heart of quantum feedback control. 

Theorem 4. The expectation J^. j[ ) of the operator valued cost when a feedback 
control strategy x r is in operation can be written as a classical expectation 

(4.7) J r ( ?) x) = / J?( c , x)P (dw) = En[j;(?, x)] 



0/ t/ie random cost-to-go function 

(4.8) j?( C ,x) = / T (CM| C (^^]))dt+(?^ r HI«) 

where q l ~^ (w) ?s f/ie solution q (i, w) to i/ie controlled filtering equation correspond- 
ing to the chosen measurement process Y\ classically represented as for the 
feedback strategy k with the initial condition q (r) = 

Proof. Using the existence and state invariance of the conditional expectation and 
the classical isomorphism proved in the first Section, it is straight forward applica- 
tion of the formula (|4.6|) to x (t) — c (u (i)) and 2: (T) = s in 14. II We can write the 
expected cost as 

Jr(?,Xr) = / (?®e T -|7r,« r (c(w(i)))>clt+(?®e r |7n« r (s)) 

T (tu [r , e* [7^ (c(« (t)))] ) di + (tn [r |e^ [ 7 ^ r (s)] ) 
T E [(e^|c( M (i)))]dt + Ej(^ u n^ 



□ 

This allows us to treat the quantum control problem as a classical control problem 
on the space of quantum states. We define the expected cost-to-go by the classical 
expression (|4.7[) for a truncated admissible strategy u r S /C r when starting in 
an arbitrary state q at time r where c r ~^ (w) is evaluated at w 6 f2 4 i solution 
q (<) = <^ r ~^ (?) to the controlled filtering equation for these initial conditions at 
t = r and the initial strategy x r for r < t < T. 

Theorem 5. Suppose that S(t, ?) is a functional which is continuously differentiable 
in t, has continuous Frechet derivatives of all order with respect to q and satisfies 

(4.9) u mf yt) {(,,c( M )) + EM^S(U)]}=0. 

for all < t < T and S(T,q) = (q,s) /or aH q G <S. Suppose also that x° is the 
strategy built from the control laws attaining these minima within a convex space 
U of admissible control values, then S(t, q) is the functional which minimizes (|4.7[) 
and x° is the optimal strategy for the control problen^. 



2 Additional technical assumptions and mathematical rigour are required to formalise the proof 
of this theorem when dealing with unbounded operators which is beyond the scope of this paper. 
See recommended texts e.g. |64l.|65| for a formal classical treatment. 
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Sketch proof. Let {it (i)}, {^ol be any control and state trajectories resulting from 
an admissible strategy >c on the initial state then from (|4.9|) . we have the in- 
equality 

(^,2(u(t)))+En[^S(t,?* )] >0 

which we integrate over [0, T] and take the expectation En to obtain due to con- 
vexity 

En [ / (? , c(u (t)))dt + S(T, - S(0, > 
Jo 

Since we have £q = ? initially, we can rearrange and use the terminal condition to 
obtain 

S(0, < En [ / (4 2(« (*))) d * + <?o . «>] = Jo(c, *)• 
Jo 

with equality when k = k° and so the lower bound S(0, is attained, proving 
optimality of xr°. □ 

A choice of controlled filtering equation is required to determine the stochastic 
trajectories d? along which to differentiate candidate solutions. The next two sec- 
tions are concerned with the examples of feedback control with respect to QND 
measurements of the diffusive process an( i the counting process {^jy} re- 

spectively. 

4.2. Quantum state Bellman equations. First let us introduce notations of dif- 
ferential calculus on the quantum state space 6 C b + . Let F = F [•] be a (nonlinear) 
functional <r i— » F [?] on 6 , then we say it admits a (Frechet) derivative if there 
exists an b-valued function V<;F [•] on b* such that 

(4.10) km \ {F [• + hr] - F [•]} = (r|V f F [•]) = (r, V,F [•]) 

for each r = r* e b*. In the same spirit, a Hessian V® 2 = V.; <g> V ? can be defined 
as a mapping from the functionals on to the b^j^-valued functionals, via 

hSsto ~hb {f[ ' + hT + H ' T ' ] - F [• + H - F [• + h'r'} + F [•]} 

(4.11) =(r«)r'|V s ®V f F[-]}. 

and we say that the functional is twice continuously differentiable whenever V® 2 F [•] 
exists and is continuous in the trace norm topology. 

With the customary abuses of differential notation, we have for instance 

v ? /« ? |x» = /'«<r|x»x 

for any differentiable function / of the scalar x = (s\X). 

4.2.1. Diffusive Bellman equation. We have the Ito rule dF/dY"^ = Sikldt with 
i, k G I e = Iw for the increments of the diffusive processes Y* = v t (W*) which 
have the expectations E^kdY"/] = (q\2$t(Li))dt, so using the classical Ito formula 
for the diffusive process Y? we can show 

Ej,[As(t,0] =^S(t, ? ) + (<r|A(«)[V t S(t,0]> 
(4-12) ' i ^ 

+ 2 E (^(0®^(0|Vf 2 S(t,0) 
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where V<;S G b denotes Frechet derivative with respect to ? £ b*, and V® 2 S G 
b® 2 := b<g>b denotes Hessian applied to S. Observing J^S(i,?) is independent of u 
leads to the following corollary. 

Corollary 2. Suppose there exists a functional S(t,s) which is continuously differ- 
entiable in t, has continuous first and second order Frechet derivatives with respect 
to ? and satisfies the following Bellman equation 

(4.13) -^ S ^) = l%{mn) + X(u)[W,S(t,,)})} 



J $>'(?) ®<T(?)|Vf S(i,0> 



2 

for all t > 0, s E S with the terminal condition S(T, ?) = (?, s). Then the strategy 
k° (^t, Y'^j = u° (i, ?) built from the control laws 

(4.14) u°(t,0 =arginf \{<;\c(u) + \{u)[V,S(t,<;)]}} 

for < t < T is optimal for the feedback control problem based on diffusive output 
measurements. 

Note that last line in 14.131 is precisely half of the Laplace operator 

asm= £ (J l (0®<rwivf 2 s(^)) 

in the quantum state 'coordinates' ? G b* as the sufficient coordinates from the 
preadjoint space of the algebra b. 

4.2.2. Counting Bellman equation. We have the ltd rules dY*dY£ — 5l5 J k <lYj for 
the increments of the counting processes Yj — v (iVj), j G I e '■ In which have the 
expectations E^[dY"/] = (s\Li (t) L* (t))dt, so using Ito formula for the counting 
processes K* we we can show 



E( 2 



=^S(t,0 + <c|A(u) [V,S(t,0]) 



( 4.i5) - E fci L * L *> (« l w - ? l V < S ^ o> 

+ £ <^*}(S(i,a l (?))-S(i,?)) 

The last two lines can be written in the Feller form as ^/u, (?) A 1 S (t, ?) as in the 
diffusive case in terms of doubled difference combination [i { (?) A 1 S (t, ?) of the 
linear combination of differences 

S (*,«*(?)) -S(t,?). Mi (<0 = <<^*> 

and //j (?) (5* (?) |V ? S (£,?)) in terms of <f (?) = a 1 (?) — ?. The formal Taylor 
expansion 



oo 

A l S(t,?) = 2^-(<f (?f"|Vf"S(t,?) 



n 

n=2 



of each A l S (t, ?) = 2 (S (a (?)) - S (?) - (<5 4 (?) |V S )) in terms of the higher order 
Frechet derivatives V® n starts from the Hessian V® 2 = V ? <S> V ? , determining the 
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Laplace operators ^<T (?) \Vf ) for the diffusive approximation of this counting 
measurement case. 

Let us introduce the Pontryagin's 'Hamiltonian' in the 'coordinates' q (?) = g — <r 
and 'momenta' p € b as the Legendre-Fenchel transform 

H (q,p) = sup {(\(u) [q]\p)-L(q 7 u)} 

uev 

of the 'Lagrangian' L (q, u) = (g — q\c (it)), where g = q—s is any stationary element 
g = g* in b* such that A (u) [g] = for any u E U (e.g. g = 0). Then one can 
write Bellman equation defining minimal expected cost-to-go (|4.5[) as the action 
functional in the compact Jacobi-Feller form as follows, similar as it was done for 
the diffusive case. 

Corollary 3. Suppose there exists a functional S(t,<;) which is continuously dif- 
ferentiable in t, has continuous first order Frechet derivatives with respect to ? and 
satisfies the following Bellman- Jacobi-Feller equation 

(4.16) -^S(t,c) + H($(0,V,S(t,0) 

(4-17) = i ]>>|L 4 L*)A*S(i,0 

for all t > 0, q € b* with the terminal condition S(T, <;) — (<;, s). Then the strategy 
h° (t, ^Tv^j = u° (t, <;) built from the control laws 

(4.18) u°(i,c)=arginf {(^(«) +A(«)[V,S(U)])1 

is optimal for the feedback control problem based on counting output measurements. 

Thus we have shown that without loss in optimality, one can reformulate the 
unobservable quantum feedback control problem into a feedback problem based on 
indirect QND measurements with feedback of the controlled conditional density 
matrix. However, the corresponding Hamilton- Jacobi-Bellman equation resulting 
from the minimization rarely has a regular solution S(t, c,) from which to construct 
the optimal feedback laws. We now study a specific quantum filtering and feedback 
case where such a control solution can be explicitly found, which is familiar as the 
only such example in the classical case. 

5. Application to a Linear Quantum Dynamical System 

Let x m — (xi,...,x m ) be the row-vector of self-adjoint operators Xj = x*, 
j = 1, ...,m and J = (Jifc) be an anti- symmetric real valued matrix defining the 
canonical commutation relations (CCRs) 

(5.1) [in, Xk] ■= Xiik - x l ii — ifiJifcl, 

written in the matrix form as = iftJl where 1 is the identity operator 

on f). Usually this is the standard symplectic matrix J T = —J = J -1 , but we 
may not assume that J is standard or non-degenerate in order to include also the 
commuting random variables x m = (xi, . . . , x m ) as a special (classical) case. It is 
worth remarking at this point that the noncommuting operators Xj are secondary 
commuting (in the sense of commutativity with all the commutants (|5.1| ). and 
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therefore they must be unbounded in the Hilbert space fj, affiliated to the the 
generated algebra b. 

We couple the open quantum system to 4 = \I e \ estimation (side) channels 
with linear combinations Li = y\ AgXj indexed by a subset I e C {1, ...,d}, 
given by a complex- valued d e x m matrix A e = (Ay) with A, y - = for i ^ I e . 
In the i-th estimation channel, we perform a measurement of the output Yf — 
(23? [Ar (t)]) defined by v t \A = v^A^ where 23? [At] = Af + A l _, i g I e 
and the output transformations is given by a quantum stochastic evolution v t on 
the algebra b (g> Ao generated by the canonical independent variables x, and Af . 
The system is also coupled to df — feedback (input) channels by the operators 
K l = L* as linear combinations of (xj), given by the row Lf = i.Aj of operators 
Li, i g If for a complex- valued df x m matrix A/ = (A y ) with Ay = for i ^ If, 
and we apply input controls with real-valued components Ui (t) 6 1 in the row 
u = (ui) ieIf via the i-th feedback channel at time t. Both matrices A e , Af may 

depend on t, but they are always orthogonal such that A^A = AgAe + A^A/, where 
A = A e + Af and A^ is transposed to complex conjugated matrix A* = (Ay) . 

Let the free dynamics of the quantum system be described by a quadratic Hamil- 
tonian Ho — ^i;.M _1 i;I for a symmetric real m x m matrix M _1 . We now introduce 
a coherent control source separating the controls u(t)Hfdt in coherent superpo- 
sition with the noise S3 (d^4+) B/ in control channel, where 2i3 (Af) — Af — A l _, 
coming from the feedback channel with B/ = 2 Re Ay, so that the Hamiltonian in 
(|4.4j) is modelled by 

(5.2) H 2 u(t) = ~x.M- 1 xZ+u{t)B f xZ. 

Using CCR's (|5.1I) and assuming that the QS dynamics is purely diffusive (no 
scattering, S l k — 18\ in (|3.13[) ). we can easily evaluate Lindblad generator A„ (xf) 
for controlled CP hemigroup dynamics f r u (i, x) — t\_ t (it, x) arising from time 
dependent quadratic Hamiltonian (|5.2p . Substituting Xj into the decomposed gen- 
erator (|4.4[) with linear L l in Xj and zero jump part, p£ (x) = 0, we obtain \ u (xf) 
as linear transformation of x, — (£i) written in vector form as 

X u (t,xl) = J(M" 1 +hlm (A f A))^I -CfU 1 (t) 

where = By J. From this, we obtain the quantum Langevin equation for 
X, (t) = vt (£• <£> I) in the linear form 

(5.3) dX. (t) + (x. (i) A T + u (t) Cj) dt = dU l , 

derived in [55J, where A T := (ftIm(A T A*) + M -1 )J with quantum noise 

(5.4) Ul := 2H% [Af (t) A] J = + W.*, 

for A = A e + A/. Here Vf = VfC T is total Langevin force and W* = W*FJ is total 
Wiener noise given by left action of matrices 

(5.5) C T = 2ReAJ, F T = fi.ImAJ, 

on rows V l = 23 (Af , ...Af) (t), W* = 23? (Af, . . . Af) (t) of all quantum Langevin 
forces and all conjugate Wiener noises respectively, coming both from the estima- 
tion and feedback channels. 
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5.1. Quantum filtering of linear, Gaussian dynamics. The linear output 
equation for the row Y e — {Y*) ieI of observable processes in estimation chan- 
nel satisfies 

(5.6) AY\ = X. (t) B T e dt + dW* e 

where B e = 2 Re A e . The quantum measurement noise is given by the row W e = 
(Wi)i£i having Gaussian independent increments on each measurement channel 
with zero mean and standard variance given by d e x d e identity matrix I e . Con- 
sidered alone, this noise represents the standard c? e -dimensional classical Wiener 
process which we measure in the field after interaction with the quantum object by 
the output isomorphic transformation Y* = 25ftuoo [A* (<)] for = . 
However, it does not commute with the quantum Langevin force V* as 

(5.7) lV;\W s e } = (rAs)ihCJ, 

due to the noncommutativity with V\ = (V*) ieJ and commutativity with = 
H$s A'j (t) , resulting form independence of A~j = (Afj.j . The fundamental 

CCR (|5.7p . defined by nonzero (if J ^ 0) matrix CJ = B e J, was first derived in 
a complex form in 18J, 44J, and in even more general infinite dimensional setting 
in |28j . It expresses the Heisenberg error-perturbation uncertainty principle in a 
precise form 

h 2 

(5.8) dVJdV. > —C e Cldt, dW e dWl = I e dt 

derived by Belavkin in [44j . [28j as necessary and sufficient condition for nondc- 
molition causality of the observable past Y l \ = {Y r e : r < t} and quantum future 
described by {X, (s) : s > t}. 

Usually in classical filtering theory the process U* and measurement noises Wl 
are considered to be independent, although in the quantum setting the Heisenberg 
principle, implying the dependence of XJ\ = V* + Wl and W*, may result in a 
nonzero covariance matrix F e = SJ Im (A|j describing the real part of quantum 
ltd table 

(5.9) dUJdW e = (f e + jC e ^j dt 

as the sum of dV. T dW e = f C e dt and dW?dW e = F e dt. Note that although each 



component Ui = U* of the row U, = (Ui, . . . , U m ) for vector quantum noise (|5.4p 
having the independent increments can also be realized as a classical Wiener process, 
these components mutually do not commute, having complex multiplication table 
dU]dU, = f^A^Adt with imaginary part defining the commutation relations 

[C/; T , U^} = (r A s) 2i/i 2 J T Im (A f A) J. 

The symmetrized multiplication 5ft [dt/JdC/.] = h 2 Re [At A] dt results in the sym- 
metric covariance 

Re (V; J V. S ) = (r A s) h 2 3 T Re (At A) J T 

defined by At A = A|A e +AyA/. It can be parametrized as FJF e + G with positive 
matrix 

(5.10) G= ^.CeCJ + ^JTRe^AjA/) J 
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implying the error-perturbation uncertainty relation 

SR [dUJdU.] > (pjF e + (h/2) 2 C e Cj) dt 

in terms of the total perturbative noise (|5.4| in the Langevin equation (|5.3p with 
respect to the standard normalized error noise in estimation channel (|5.6|) . 

Let us denote the initial mean (x.) of the phase space operator vector by the 
component wise expectation x, = (s,x») and covariance matrix 

S := (Re(<r, XiX k ) - xjx k ) 

where Re(<r, XiXk) = 5^, ijXfc + XkXi), which is symmetric positive definite matrix 
S = (Ejfc) satisfying the Heisenberg uncertainty inequality 

(5.11) £>±-J 

As it was shown by Belavkin in {5j , |66j , and even in infinite dimensions in [28j , 
the filtering equation ()3.24|) preserves the Gaussian nature of the posterior state 
[28j . so the posterior mean x\ = (<: ,x a ) and the matrix S (t) of symmetric error 
covariances as real part of 

(?o> ~ (^fe ~~ %)) = (€o,%i%k) — x\x\ 

form a set of sufficient coordinates for the quantum LQG system which agree with 
the initial mean and covariance for q (0) = <;. Applying (|3.24|) to the first and second 
symmetrized moments of x m , i.e., rigorously speaking, to the spectral projectors of 
these unbounded operators affiliated to b = B (fj), provides posterior expectations of 
these sufficient coordinates for diffusive non-demolition measurement of the output 
operators Y l e . These can be found as solutions to Belavkin's Kalman filter equation 
[44], [28] written in vector form as 



(5.12) di\ + (x^A 1 + u (t) Cj)dt = dW^K T (i) 
with initial condition x® = x, for the posterior mean, 

(5.13) K(i) = £(<)BT + F e , dw[ = dYl - x*.Bldt 

and for the symmetric error covariance we have 

(5.14) -^E = G - SAJ - A e S - SBm e S, 
dt 

(5.15) A e = A + F e B e , S(0) = S. 



5.2. Quantum LQG control and microduality. We aim to control an output 
quantum stochastic linear evolution Z (t) = v Q {z ® I) of a d /-dimensional linear 
combination z — i.Ep where E/ is a real, in general time dependent df xm matrix, 
represented in the Heisenberg picture as a row [Zj) i&1 — A.Ej of Zi (t) = v t (i,) by 
forcing Z (t) to follow the classical input trajectory of the feedback controlling force 
u (t) whilst constraining for energy considerations a positive quadratic functional 
of phase space operators X, (t) = vt (x, ® /). Thus, our control objectives and 
restraints can be described by the general quadratic operator valued risk (14. ip in 
the canonical form 

(5.16) c(u) = (u - z) (u - z) T + x.Hx T , 

and s = x t £lx~l for positive real symmetric m x m matrices f2, H. 
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Since x, and 53 are generators for the full probability distribution given by the 
Gaussian posterior state they form a set of sufficient coordinates, so we may 
consider derivations of S(t, as partial derivative of S(t, x a , 53) 

(<k, V c S(t, 0> dx.dJSit, x., S) + (d53, d"S(t, a:., 53)) 

+ ~(ft®ft,Vf 2 S(t,0> ~ +i(da;Tda;.,aTa - s(t,a:. ) S)). 

We use the notation (•, •) to denote the matrix trace inner product (D,F) = 
Tr[D T F] on the vector space of matrix configurations for the multi-dimensional 
system. This gives the directional derivatives along dxj and dTiik as functionals of 
the column djS of partial derivatives d l S = dS/dxt and the matrices 

d S (<953 ii; S ) ' ^.d.S i^dxidxk^ 

which are evaluated at (a:., 53 (i)). 

Inserting this parametrization into the Bellman equation (|4.13p and minimizing 
gives the optimal control strategy 

(5.17) u(f) = -d'S(t,:r*,53)C+:E*.Ej 

where S(i, x.,53) at x, — £^,53 = 53 (t) now satisfies the non-linear partial differ- 
ential equation 

-^S(t,ar.,E) = l(x,A^dlS + d.SAxl) + x.Hxl 

+ (G-AS + SAT,a"S) + (5:,H) 
l& 8j -(|a.SC + x.E})(ia.SC + x.E})t 

+ ((53B T + F e )(53B T + F e ) T , ffid.S - d"S) 

which is the Hamilton- Jacobi-Bellman equation for this example. 

It is well known from classical control theory that LQG control has a minimum 
cost-to-go which is quadratic in the state, so we try the candidate solution 

S(t, x„ 53) = x.n(t)xl + («(*),£) + a{t) 

in the HJB equation (|5.18[) . This separates the HJB equation into a set of coupled 
ordinary differential equations and gives the optimal feedback control strategy 

(5.19) u(t) = x. (i)L T (f), L T (t) = tt(t)Cf + Ej- 

which is linear in the solution to the filtering equation cc* at time t where ft(t) 
satisfies the matrix Riccati equation 

(5.20) ~T° = H ~ OA / -AJfi-fiC/Cjft 



dt 

A f = A + C/E/, O(T) = $7 



and a(t) satisfies 
d 



(5.21) --a(t) = ((n(t)C / + E})(n(t)C / +E})T,S(t)) + (n(t),G) 
a(T) = 0. 
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From this we obtain the total minimal cost for the control experiment 

(5.22) S(0,a;.,E) = x.floxl + Tr[« £] + / Tr[ft(i)G]d< 

Jo 

+ [ Tr[(ft(i)C / +Ej) T £(i)(ft(t)C / +Ej)]ctt 
Jo 

where f2o is the solution to (|5.20|) at time t = 0. 

The equations (|5. 13[) - (|5. 14|) and (|5 . 19|) - (|5 . 20[) demonstrate the intrinsic duality 
between optimal quantum linear filtering and optimal classical linear control, which 
we call microduality. To make this duality more transparent let us introduce real 
matrices E and F defining the matrices F e and E/ in (|5.5p and (|5.16|) by EJ = FJ 
and E/J = F T in the similar to BJ = CJ and B fj — C T in terms of the estimation 
and feedback channel matrices B = B e and C = C/. Then the microduality is 
summarized in the table 



(5.23) 



Filt 


AJ 


BJ 


EJ 


K (T — i) 


G(T-t) 


£ (T - 1) 


Con 


JA T 


C T 


FT 


JL T (t) 


JH(t) J T 


jrj (t) j t 



in which the duality notations are made in filtering-control alphabetical order 
(B,C), (E,F), (G,H) and (K, L) and the matrices A B,E should be also taken 
at t T = T — t if they depend on t for the duality with A T ,C T ,F T evaluated at 
t. This duality allows us to formulate and solve the dual classical control problem 
given the solution to quantum filtering problem with dual parameters. The duality 
can be understood when we examine the nature of each of the methods used. Both 
methods involve the minimization of a quadratic function for linear, Gaussian sys- 
tems, (i.e. the least squares error for filtering and the quadratic cost for control). 
The time reversal in the dual picture is explained by the interchange of the input 
(feedback) and the output (estimation) channels together with the linear canonical 
transformation given by the symplectic matrix J. Note that E = fiImA e , and G 
must be positive definite satisfying the relation (|5.10|) due to Heisenberg uncertainty 
relation corresponding to the error-perturbation CCR's [dWe,dZ7 e ] = iftF e CJdt 
in the Ito multiplication table (??). In order to complete this filtering-control mi- 
croduality we may set F T = Mm A/ J to have the relation between E and F similar 
to the duality of B = 2ReA e and C T = 2ReAJ-J, and also assume that matrix H 
is also positive, satisfying 

(5.24) H = (ft/2) 2 BJB/ + h 2 Re (A^A e ) , 

where B/ = 2 Ira Af. Note that although the condition (|5.24[) is not a requirement 
in this classical-quantum setting, in which only positivity of the combination H + 
EjE/ suffices, this relation may be required for the fully quantum setting when the 
both the input it and the output Y e are allowed to be noncommutative, which will 
be studied and published els where. 

6. Discussion 

The Bellman equations for quantum systems having separate diffusive and count- 
ing measurement schemes have been derived in continuous time under a general 
setup. This presents original derivations of the results stated in [5] (a derivation 
of the diffusive case has also appeared in [43]) and allows us to reformulate the 
optimal control problem for a fundamentally unobservable quantum system into a 
classical control problem on the Banach space of observable filtered states. 
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The multi-dimensional quantum LQG problem which finishes the paper demon- 
strates the first application of the general quantum Bellman equation from which 
one can obtain the special cases of the Gaussian quantum oscillator |5] and quan- 
tum free particle [67] HI] . Note that the fundamental difference between this ex- 
ample and the corresponding well studied classical case is in the observability of 
the system. The quantum noises introduced act only to account for the quantum 
backaction due to the incompatibility of quantum events. No further restrictions 
on the observability or additional classical noise are introduced. As such, the cor- 
responding classical problem (when h — > 0) admits direct observations of x\ = x\ 
for a deterministic classical system and has an optimal direct feedback strategy 
u (t) — x\Jj(t) and minimal cost S(0,x.) = x,SIqxI- Also this example clearly 
demonstrates the micro duality between quantum linear filtering and classical feed- 
back control as a more elaborated duality involving also the linear symplectic 
transformation J. 

During the publication of this paper, there have appeared a number of recent 
works on quantum filtering and feedback control by Bouten et al. [68l ESI EQ] 
to which we refer the interested reader for a more detailed introduction to these 
subjects. 

7. Appendix 

7.1. A. Some definitions and facts on W*-algebras. 

(1) A complex Banach algebra A with involution a i— > a* such that ||o*ffl|| = 
||a|| 2 is called C*-algebra, and W*-algebra if it is dual to a linear subspace 
LCA* (called preadjoint of A = L* if it is closed, denoted as L = A*). 
They all can be realized as operator algebras on a complex Hilbert space 
H, and an operator W*-algebra is called von Neumann algebra if its unit 
is the identity operator / in Ti. The simplest example of W*-algebra is the 
von Neumann algebra B (7i) of all bounded operators acting in a complex 
Hilbert space H. A von Neumann algebra A is called semisimple if H has 
an orthogonal decomposition into invariant subspaces Hi in which A is 
B(7ii). Let {Qi} (or {-4i}) be a family of self-adjoint operators (operator 
algebras Ai) acting in H, e.g. orthoprojectors Q\ — Qi — Q*. The W*- 
algebra generated by this family is defined as the smallest weakly closed self- 
adjoint sub-algebra A C B (TL) containing these operators, or the spectral 
projectors of these operators if Qi are unbounded in 7i (or the algebras Ai, 
A = VAi). It is not necessarily semisimple but in the case I £ A it consists 
of all bounded operators that commute with the bounded commutant B = 
{B e B (H) : BQi = QiB Vi} = {Qi}' (or with B = rU<), i.e., it is the 
second commutant A — {Qi}" of the family {Qi}. The latter can be taken 
as the definition of the von Neumann algebra generated by the family {Qi}. 
Note that the commutant B is a von Neumann algebra such that B = A', 
called the commutant algebra of A. [71] . 

(2) A (normal) state on a von Neumann algebra A is defined as a linear ultra- 
weakly continuous functional A — » C, satisfying the positivity and normal- 
ization conditions 



(7.1) 



(Q,Q)>0, VQ^O, (q,I) = 1 
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\Q ^ signifies the nonnegative definiteness (ip\Qip) ^ Vip G TC called 
Hermitian positivity of Q] . The linear span of all normal states is isometric 
with the preadjoint space „4*. The latter is usually described as the space of 
density operators <; uniquely defined as (generalized, or affiliated) elements 
of the algebra A with respect to a standard pairing (g*,Q) := (g\Q) = 
g* (Q) of A* and A given by the mass /i (g) = (g\I) on the positive g such 
that g* — g > is state density iff (g\I) — 1. A state g is called vector state 
if {g,Q) = {ip\Q4>) = {g^\Q) (denoted g = g^) for some ip 6 H, and pure 
state if it is an extreme point of the convex set S (A) of all normal states on 
A. Every normal state is in the closed convex hull of vector states g^ with 
\\tp || = 1 but there might be no pure state in S (A). If algebra A is semifinitc 
(there exists a faithful normal semi- finite trace Q t— » trQ, then the states 
on A can be described by unit trace operators g G A (or g h A if the are 
only affiliated to A), by means of the tracial pairing {g\Q) — tr [f?*<2]. In 
the simple case A = B (Ti.) the density operator g is any nuclear positive 
operator normalized with respect to the usual trace [71] . 
Let A, B be von Neumann algebras in respective Hilbert spaces TCq and 
Hi, and let $ : B — > .4 be a linear map that transforms the operators 
B € B into operators A G .4 (called sometimes superoperator) . The map $ 
is called a transfer map if it is ultraweakly continuous, completely positive 
(CP) in the sense 

oo 
i,k=l 

(i = l,...,d e < oo), and unity-preserving: $ (ii) = I (or ^(/i) < /o)- 
The CP condition is obviously satisfied if $ is normal homomorphism (or 
W*-representation) ir : B — > A, which is defined by the additional multi- 
plicativity property ir (B*B) = it (B)* it (B). A composition (g , $ (£?)) = 
(g Q \& (B)) with any state g = g$ is a state g x on S described by the adjoint 
action of the superoperator $ on j : 

(f> |$ (£)) = (#* (ft,) \B) ,V5e6,g e A- 

A transfer map $ is called spatial if 

$ (B) = FBF* or (g Q ) = F^ F;, 

where F is a linear coisometric operator Hi — ► ?ioj -F 1 ^ 1 * = -^o (or FI\F* < 
Jo) called the propagator and = is defined as left adjoint (F*p|Q) = 
(g|FQ) with respect to the standard pairings (which is usual adjoint, F» = 
F* , in the case of tracial pairing (g\Q) = tr [p*Q]). Every transfer map is 
in the closed convex hull of spatial transfer maps, but there might be no 
extreme point in this hull. 

Let V be a measurable space, and 05 its Borel tr-algebra. A mapping II : 
dv G 03 I— > II (dv) with values II (dv) in ultraweakly continuous, completely 
positive superoperators B — > A is called a transfer measure if for any tp G 
A, B G B the C- valued function 

(n(dv)*g\B) = ( Q \n(dv)B) 

of the set dv C V is a countably additive measure normalized to unity 
for B = I. In other words, II (dv) is a CP map valued measure that is 
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er-additive in the weak (strong) operator sense and for dv = V is equal 
to some transfer-map $. In particular, H(dv,B) = M(dv)$(B) with 
M (dv) = II (dv, I) is transfer map iff [M (dv) , $ (B)] = for all dv G 03 
and B G B as it is the case of the nondemolition measurements given by 
representations M = E of 03 in A and <E> = 7r in E (03)' n A. The quantum 
state transformations g i— » II* (A, g) corresponding to the results v G A of 
an ideal measurement are described by transfer-operator measures of the 
form 



where F (v) denote linear operators Hi — ► Ho, the integral with respect to 
a positive Borel measure A on V is interpreted in strong operator topology, 
and / F (v) F (v) dv = Iq. Every transfer-operator $ : B — > A can be 
represented by the integral (|7.4| on A = V of an ideal measurement II (dv) 
as the compression 



of the nondemolition measurement E (A) ir (B) = 1a <8> B on the extended 
Hilbert space H = L F (v)* Hodv with the isometric embedding F*-0 O = 
f® F (v)* tp X (dv) of ip G Ho into H adjoint to the coisometry Ftp = 
/ v F (is) ip (v) A (dv) of V G H into H . 

. B. Notations of quantum stochastic calculus. 

(1) Guichardet Fock space = T(£ t s ) over the space £ t s = L 2 (If) of square 
integrable complex functions £ (r) on the interval If = [t, t + s) is built as 
the Hilbert sum ©^L £n °f the spaces £„ = L 2 (r„) of square integrable 
functions x n : t G T„ i— > C of finite chains r = {ti < . . . < t n } G If 
for all n = 0, 1, . . . where Xo (0) = c is a constant corresponding to only 
one - empty chain To = of To, with c = 1 for the vacuum vector state 
Xn = 5q ■ It is L 2 (r (If)), where T (If) is disjoint union £^° =Q T„ (If) of the 
n-simplices T„ (If). The integration on T is assumed over the Lebesgues 
sum dr = dr n of measures dr„ = dti . . . dt n on the simplices T n of 
chains r„ : r„| = n with the only atom dr = 1 at T = {0} such that 



for the exponential vector-functions £ (r„) = £ (ii) • • •£ (*n) given by a 
single-point function £ G £. It is isomorphic to both usual Fock spaces 
©nLo^ °^ symmetric and antisymmetric functions x„ ( r i, ■ ■ ■ , r n) extend- 
ing x ( r «) on (If)™ with respect to the measure (nl) dri---dr„. Fock 
space is infinitely divisible in the multiplicative sense T\_* ~ Tl_ r ®T^ for 
any r, s > 0, which is a reflection of the additive divisibility £l_? ~ £f_ r ©£ t s 
of L 2 ) • The generalization to the multiple Fock-Guichardet case over 
the space £* — L 2 (If — > {!) of vector- valued functions £ (t) in a Hilbert space 
£ (= C d , say) is straight forward and can be found in [54], [24]. All prop- 
erties remain the same, and the only difference is that T\ is not L 2 -space 
of scalar-valued functions \ on T but is Hilbert integral J^ r ^ s ^ IC (t) dr 




H(A,B) = FE(A)tt(B) F* 
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of K.(t) ~ £®I T I ®L 2 (I(), the spaces of square integrable tensor-valued 
functions x : T„ -> 6® n . 
(2) Four basic integrators At, A°_, and A° of the universal quantum sto- 
chastic (QS) calculus [54], [24] are operator- valued measures A v ^ (I) of preser- 
vation, annihilation, creation and exchange respectively, defining the basic 
QS integrals of the total QS integral as sum-integral 

(7-6) i* (*0= E f*KS(r)A^dr) 

|U=-,o;i/=o,-|- *' 

of four basic integrants if 7 , K~ , if? and if ° as measurable operator- valued 
functions if£f (r) in J 7 by the following explicit formulas: 



(7.7) 

(7.8) 
(7.9) 

(7.10) 



(7.11) 



(7.12) 



X(r) = 



\il(K)°_]x(r) 
%(K)tx] (r) 

[i(KTo]x(r) 



[if [r) X ] (r)dr, 



[if (r) X (r)] (r) dr, 
]T^Mx] (r\r), 



Here [% (r)] (r) = x (?" LJ t), where rUr is union of disjoint r 6 T and r ^ r, 
r\r is difference of r and a singleton r = {r} C r and Tg = r n Ig. The 
functions if£f (r) should be L p -integrable in a uniform operator topology 
[51], [21], with p = 2j[v-n) where - = -1, o = 0, + = 1. Note that 
these definitions do not assume adaptedness of integrants as they generalize 
Hitsuda-Skorochod extended stochastic integral. The multiple version of 
this explicit QS-integration is straight forward and can be found also in 
[54] , [24] , and the adapted version based on coherent vectors is in [50] . 
(3) Ito rule (O for QS integrals M (t) = M (0) + i& (if) with adapted four- 
integrant K (t) = (if£) (t) is based on the noncommutative Ito table (|3.10j) 
which uses ^-matrix algebra of the canonical triangular representation 



K = 



K- K- 
K° if ° 




if* 



R- R+ 
R° R° + 




where R t t v = K'L*, for any noncommutative ltd algebra [54] . [24] . It was 
derived in 48J for simple bounded integrants, and therefore can not be 
rigorously applied for multiple integration of quantum stochastic equations 
except the special unitary case. In the general form presented here QS 
Ito formula was proved for unbounded integrants in [53] where it was also 
extended to nonadaptcd integrants, and the functional noncommutative ltd 
formula was also obtained in the pseudo-Poisson form as 

df (M (t)) = (/ (M (*+)) - / (M (*)) ® I f v A; (dt) , 

where M (t+) = M (t)®I + K (t) is QS germ [72] of the QS integral M (t) 
which is defined by its four QS-derivatives Ktf (t) and unite matrix I = (5„), 
and the summation convention over /i, v — — ,o,+ is applied. Using this 
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formula the HP differential conditions (|3.12j) of QS unitarity of QS interac- 
tion evolution Ut were obtained as pseudo-unitarity condition in terms of 
germ U t+ — Ut® S, and also QS differential conditions of complete posi- 
tivity, contractivity and projectivity were found in |72j,}58] respectively as 
its pseudo complete positivity, pseudo contractivity and pseudoprojectivity 
of the corresponding QS germs M (t+). 
(4) Using quantum Ito formula the general QS evolution equation for a quan- 
tum stochastic density operator g (t) was obtained in the form of quantum 
stochastic Master equation [55]. [50] 

(7.13) dg (t) = (Gig (t) G? - g (t) S L K ) dA?, g(0) = g 

Here G = I + L is the germ of QS evolution equation dVt = VtL^dA^ which 
is not assumed to be pseudo-unitary. In the case of normalization condition 
S~ = O in terms of left adjoint operators St v = G_„ with respect to a 
standard pairing (b*|b), this equation is called QS decoherence, or entan- 
gling equation for quantum states satisfying normalization {(q it) , I)) % = 
(g,i) with respect to the pairing (g*,A)^ — (gS^\A6^) = (glA)^ on the 
noise algebra A given by the vacuum vector 5® = 0®. This is the most 
general QS equation preserving complete positivity and normalization in 
this mean form. Denoting K~ = —G~ = K L such that G*+ = —K*, this 
can be written [59], [58] as 



dg (t) + 25R [K L g (t) dAt_] = ]T G 3 K g (t) G{* - g (t) ^ dA« L 



More explicitly this Belavkin equation can be written in terms of K = K+, 
L\ = G\ = L l such that G 3 * = V* as 



dg (t) + 25R [Kg (t)] - V Q (*) V * dt 



5> 



Y,G{g{t)Lr* -K k g{t)\dA k L 



E 



Y,G{g(t)Gr -g(t)Si)dAH 



The weak normalization condition can be written as L + L* + ^2n L>iL* = 
in terms of left adjoint L.Li to —K,L l such that K = —L^,L l = L\ 
(K = —L*,L l = L* in the case of trace pairing) for any number of i's, 
and arbitrary Ki, G\, i,k = 1, . . . , d. This is QS generalization of Lindblad 
equation [57] given by the generator (|3.16[) corresponding to the case d = 0. 
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